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Through second order in perturbative general relativity, a small compact object in an external 
vacuum spacetime obeys a generalized equivalence principle: although it is accelerated with respect 
to the external background geometry, it is in free fall with respect to a certain effective vacuum 
geometry. However, this single principle takes very different mathematical forms, with very dif¬ 
ferent behaviors, depending on how one treats perturbed motion. Furthermore, any description of 
perturbed motion can be altered by a gauge transformation. In this paper, I clarify the relationship 
between two treatments of perturbed motion and the gauge freedom in each. I first show explicitly 
how one common treatment, called the Gralla-Wald approximation, can be derived from a sec¬ 
ond, called the self-consistent approximation. I next present a general treatment of smooth gauge 
transformations in both approximations, in which I emphasise that the approximations’ governing 
equations can be formulated in an invariant manner. All of these analyses are carried through second 
perturbative order, but the methods are general enough to go to any order. Furthermore, the tools 
I develop, and many of the results, should have broad applicability to any description of perturbed 
motion, including osculating-geodesic and two-timescale descriptions. 

PACS numbers: 04.25.-g, 02.40.-k, 04.40.-b, 04.20.-q 


I. INTRODUCTION 

In general relativity, a point test mass travels on a 
timelike geodesic of the spacetime geometry it resides in, 
satisfying 


D 2 z> i 

dr 2 


= 0, 


(1) 


where z^(t) is the particle’s coordinate position, r is 
proper time as measured in the spacetime’s metric g^„, 
if ;= u M := and is the covariant deriva¬ 

tive compatible with g^. 

Now suppose we take a step beyond the test-particle 
approximation and treat the mass as a small, extended, 
gravitating object surrounded by a vacuum region. The 
mass creates a nonlinear perturbation h^ atop the vac¬ 
uum background metric g llu . This perturbation acceler¬ 
ates the object, and it no longer moves on a geodesic of 
g^. Furthermore, it also does not move on a geodesic of 
the full spacetime’s metric g pl/ = g^ + h^. However, if 
the object is sufficiently similar to the structureless point 
mass—uncharged, slowly spinning, nearly spherical, and 
compact— -then it does travel on a geodesic of an effective 
metric In this case, Eq. ([l]) is replaced by [Tj 


£> 2 z^ 
df 2 


= <?(e 3 ), 


( 2 ) 


where e is a small quantity proportional to the object’s di¬ 
ameter d and mass to, and f, gy := and u M := 

are now defined with respect to g^. What is the effec¬ 
tive metric g ^It is a certain smooth vacuum solution 
g ^„ + v , where the regular field h^ v = ft M1/ — ft® is 

obtained from ftby subtracting a certain singular field 
ft®„, which encodes the local information about the ob¬ 
ject’s mass and multipole structure and diverges on zffi 


These results are derived in gravitational self-force the¬ 
ory, in which one examines the limit e —> 0 and constructs 
an asymptotic expansion of the full metric garound the 
background g^. The titular self-force in this theory is 
the force that accelerates the object with respect to the 
background metric; I refer the reader to Refs. [2H3] for 
the seminal work in this field and Refs. HUM!, together 
with the reviews US El, for rigorous and comprehen¬ 
sive treatments. Because Eq. ([2j) tells us that all objects, 
regardless of internal composition, fall freely in a vac¬ 
uum gravitational field, it can be loosely thought of as a 
generalized equivalence principle. 

However, this simple principle masks two important 
subtleties: how the form of the equation of motion and 
its solutions, and even perturbation theory as a whole, 
depend strongly on one’s basic treatment of perturbed 
motion; and how perturbed motion depends inherently 
on one’s choice of gauge. 


In self-force theory, various descriptions of perturbed 
motion are in use. In the Gralla-Wald formalism devel¬ 
oped in Refs. mi, one expands both ftand z M in 
powers of e. The equation of motion ([2]) then takes the 
form of a sequence of equations, given by (154]), (601, and 


(64) below; these are evolution equations for the zeroth- 
orcler worldline Zq and for deviation vectors that describe 
the object’s movement away from Zg. Because the devi¬ 
ations eventually grow large (due to dissipation of en¬ 
ergy through gravitational radiation, for example), this 
formalism is restricted to spacetime domains of finite, e- 
independent size. An alternative description of perturbed 
motion is provided by the self-consistent formalism de¬ 
veloped in Refs. m in izi, which expands but not z M 
in the limit e —> 0, such that z M retains its dependence 
on e. The equation of motion then takes the form Q 
given above, or in terms of background-metric quanti¬ 
ties, the form (42) given below. By avoiding the expan- 
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sion of z 11 , this formalism constructs an approximation 
valid on spacetime domains of asymptotically large size, 
such as ~ 1/e. Other treatments, such as the osculating- 
geodesics approximation mm and the two-tinrescale 
expansion of the Einstein equation nu, are intermedi¬ 
ate between Gralla-Wald and self-consistent approxima¬ 
tions, utilizing expansions of the worldline to build a self- 
consistent approximation. In this paper I will not discuss 
these intermediate treatments, as they are easily under¬ 
stood on the basis of the extreme two; I refer the reader 
to Ref. m for a description of how that understanding 
arises in the case of the osculating-geodesics approxima¬ 
tion. 

Now, supposing we have adopted one or the other 
treatment, consider the effect of a gauge transforma¬ 
tion. The test-particle equation ([Tj) is manifestly invari¬ 
ant under a diffeomorphism p: z^ and g transform to 
z'^ = p{z^) and g'^z') = p*g^, v (z) (where p* is the 
push forward), and the equation of motion is the same 
in terms of these new quantities as it was in terms of the 
old ones. Of course, the field equation is also invariant: if 
is a solution to the Einstein equation, then so is g 1 . 
Can we say the same about the perturbed equation of 
motion ([2]) and the field equation for the effective metric 
g^v that appears in that equation? Several issues arise 
when answering this question. A gauge transformation is 
induced by a near-identity diffeomorphism p, and so the 
transformations = </?(z M ) and g' = p*g^ v should still 
apply when expanded in the limit e —> 0. In coordinates, 
a gauge transformation reads x 11 —► + 0(e 2 ), and 

hence one expects —> z 11 — e£ M + 0(e 2 ); however, the 

form of z 11 , and therefore the effect of p on it, depends 
entirely on which representation of perturbed motion one 
uses, and this plays out in nonobvious ways in the metric. 
Furthermore, a transformation law for g^„ does not imply 
a transformation law for the effective metric g /il/ . How 
does this particular piece of the total metric transform? 
That is, if h is split into h and in a particular 
gauge, how does that split behave under a gauge trans¬ 
formation? 

In this paper, I present a unified framework to tackle 
these questions. I first clarify the relationship between 
the self-consistent and Gralla-Wald formalisms. In par¬ 
ticular, I show explicitly how the latter can be derived 
from the former. I then describe how, in both formalisms, 
with appropriate transformation laws for /i®„ and /ij/,, 
the governing field equations and equation of motion can 
be made manifestly gauge-invariant, just as they are for 
a test particle in a background. 

This work extends previous results in several ways. 
Primarily, it carries explicit calculations to second or¬ 
der in perturbation theory, the order at which sig¬ 
nificant complexity arises. More broadly, it provides 
very general descriptions of the self-consistent formal¬ 
ism, which has previously been formulated almost ex¬ 
clusively in the Lorenz gauge. More narrowly, in clar¬ 
ifying the relationship between the two formalisms, it 
makes concrete many ideas that had previously only 


been suggested [B] or roughly sketched PS: and in the 
case of the Gralla-Wald formalism, it obtains a covari¬ 
ant and reparametrization-invariant second-order equa¬ 
tion of motion, which had previously been obtained only 
in locally inertial coordinates 0 or in parametrization- 
dependent form m On gauge transformations, it goes 
beyond earlier work 130031HBH21 by providing the first 
clear and complete explication of the gauge freedom in 
the self-consistent formalism, for which only sketches HB3 
or incomplete descriptions m are available in the litera¬ 
ture; because this formalism is a nontrivial alteration of 
standard perturbation theory, determining the action of 
gauge transformations within it requires care. The anal¬ 
ysis also yields, for the first time, a second-order trans¬ 
formation law for the self-force, generalizing the standard 
Barack-Ori result El from first order. 

One advantage of clearly formulating gauge freedom in 
the self-consistent formalism is that it illuminates what is 
and is not invariant in long-term evolutions. The central 
practical utility of self-force theory is in modeling binary 
systems of compact objects, particularly extreme-mass- 
ratio inspirals (EMRIs) [22) [23]. EMRIs evolve on a long 
timescale of size ~ 1/e, inversely proportional to the rate 
of energy emission in gravitational waves, and accurate 
models must track the waveform on that timescale. How¬ 
ever, in formalisms designed for long-term evolution, the 
waveform, which one expects to be invariant, is a func¬ 
tional of a gauge-dependent worldline. In this paper I 
perform a preliminary analysis of this issue. Based on 
informal error estimates, I argue that while the self-force 
(and more generally, any small perturbations of motion) 
are pure gauge on domains of size ~ e°, they contain in¬ 
variant content in any well-behaved gauge on the large 
domains of size ~ 1/e, and that the gauge-dependent 
aspect of the waveform represents a small, negligible cor¬ 
rection to its gauge-invariant content. 

I organize these analyses as follows. Section [TT] sum¬ 
marizes the self-consistent and Gralla-Wald formalisms. 
Section 1III1 shows how the latter can be derived from the 
former; the methods I use can also be applied in any other 
formalism in which the accelerated motion is (implicitly 
or explicitly) expanded around a nearby worldline, par¬ 
ticularly osculating-geodesics and two-timescale approxi¬ 
mations. Section |TV| describes the gauge freedom in each 
formalism, deriving transformation laws for the various 
representations of the worldline, the physical perturba¬ 
tions h^, the singular and regular fields and h^ v , 
and the self-force. For simplicity, I restrict my atten¬ 
tion to smooth transformations. The paper concludes in 
Sec. [V] with a discussion of practical issues, particularly 
the issue of long-term dynamics. 

Before proceeding to that material, I note one subtlety 
I gloss over: the specific choice of regular held h^ v . The 
split of into /i®j, and h^ v is far from unique, and any 
number of choices can be made which preserve the core 
properties that (i) is a geodesic of g^ = g^ u + h^ v , 
and (ii), g is a smooth vacuum metric pTj [23] ■ Fur¬ 
thermore, one can also choose practical splits into h s 
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and v for which one or the other of those core prop¬ 
erties is not satisfied. For example, with the split de¬ 
fined in Ref. [ 8 ], property (i) is not satisfied, and with 
the one defined in Ref. 0, property (ii) is not. Here I 
assume the effective metric is chosen such that both (i) 
and (ii) are satisfied, but I place no further constraints 
upon it. The obvious choice is thus the one presented in 
Refs. Q] 0 EH] (a generalization of the Detweiler-Whiting 
regular field 0 from first order), but other choices could 
also be made. 

I work in geometrized units with G = c = 1 and the 
metric signature (—, +, +, +). Indices are raised and low¬ 
ered with the background metric g and a semicolon de¬ 
notes the covariant derivative compatible with that met¬ 
ric. Functionals P of a function / (or tensor field /) are 
written as either P(x; f) or P[f](x), as convenient, where 
is the coordinate label of a point. Quantities with a 
breve over them, such as z% and refer to a Gralla- 
Wald expansion; corresponding quantities without breves 
refer to a self-consistent expansion. The expansion pa¬ 
rameter e is treated as a purely formal one, to be set 
equal to 1 at the end of a calculation. Throughout the 
paper, I specialize to objects which, through order e 3 , are 
nonspinning and spherical. 

II. SELF-CONSISTENT AND GRALLA-WALD 
APPROXIMATIONS 

In this section, I briefly describe the two treatments 
of perturbed motion. This sets the stage for subsequent 
sections. It also serves to generalize past descriptions of 
the self-consistent approximation, making the formalism 
more amenable to gauge transformations. 

A. Self-consistent 

The self-consistent approximation treats the met¬ 
ric perturbation as a functional of the object’s self- 
accelerated, center-of-mass worldline. It is encapsulated 
in an asymptotic series of the foru 0 

hp U (x, e) = E z )> ( 3 ) 

7l>0 

together with an equation of motion of the form 

= F»(t, e) = £ e n F^[h 1 , ..., ft"], (4) 

n> 0 


1 All series in this paper are to be thought of as asymptotic rather 

than convergent, meaning I actually assume that for each in¬ 
teger N > 0, the perturbation can be written as h^ u (x,e) = 
^2n=o z ) + °( eN )i where “o(/(e))” means “goes to zero 

faster than /(e) in the limit e —> 0”. Generically, lne terms also 
appear in the series, but these terms do not disrupt the series’ 
well-orderedness [7]. Here I absorb them into the coefficients 
h"„(x;z). 


where are the coordinates on the worldline 7 repre¬ 
senting the object’s center of mass [a nu. The perturba¬ 
tions ft" are functionals of z ,L , and z M in turn depends 
on through Eq. (|4j). The worldline 7 is accelerated 
with respect to the background spacetime, and accord¬ 
ing to Eq. Q, its acceleration is e dependent. Hence, 
through its dependence on z^, each term z) like¬ 

wise depends on e. This e dependence distinguishes the 
approximation from an ordinary Taylor expansion. 

The idea of a self-consistent approximation has been 
around since the inception of gravitational self-force the¬ 
ory m ini, and it is by now fairly well developed D13 
0 EEU EE3 HU- However, aside from a generalization in 
Ref. [7j, the approximation has always been formulated 
in a specific gauge, the Lorenz gauge. In that formu¬ 
lation, inspired by post-Minkowski theory [23] |2B], one 
imposes the Lorenz gauge condition to split the exact, 
fully nonlinear Einstein equation into a weakly nonlinear 
wave equation plus an elliptic constraint equation. The 
wave equation is solved perturbatively for the functionals 
h™„(x;z) while leaving 7 unspecified, reducing the con¬ 
straint equation to evolution equations for the object’s 
matter degrees of freedom, principal among them the 
equation of motion governing 7 . This formulation makes 
each functional h™ v (x m ,z) a solution to a wave equation, 
with the wave operator being exactly the same at each 
order. The more general formulation in Ref. 0 allows 
one to choose a different wave operator, but even in this 
more general formulation, the wave operator is always 
the same at each order; there is no flexibility to adopt 
a gauge in which the operator takes different forms at 
different orders. 

This limitation runs counter to the usual usage of 
gauge freedom, in which one can impose a different gauge 
condition at each order and thence work with different 
field equations at each order. So we require a more 
general formulation, one which allows more flexibility in 
one’s choice of gauge. I present such a formulation here. 


1. General formulation 

To begin, first consider the exact Einstein equation 
Rfj, v [ g] = 0 in a vacuum region outside the object. After 
substituting the expansion this equation becomes 

E h n ] = E ^[ft 1 ,..., ft"- 1 ], (5) 

n>0 n>1 

where 5R^ v \j] := j^R^[g + Aj ]| a _ 0 is the linearized 
Ricci tensor constructed from a perturbation , and 
the quantities on the right-hand side, 

1 d n 

syA-.r- 1 ] --^'Ws+E^’IL.O. (6) 

p 

are made up of strictly nonlinear combinations of their ar¬ 
guments. Explicitly, S ’ 1 = 0 and = — 6 2 R IJ , l , [ft 1 , ft 1 ], 
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where is the second-order Ricci tensor; the 5 no¬ 

tation is defined in Appendix |C| 

We wish to solve Eq. 0 subject to 

(i) some desired (e.g., retarded) boundary conditions 

(ii) a “matching condition” which states that at small 
distances r from 7 , the solution takes the generic 
form of a metric outside a compact object, scaling 
as h]‘ lv ~ l/r-0 

(iii) a center-of-mass condition that makes the met¬ 
ric effectively mass-centered on 7 , as discussed in 
Refs. dJEUZUII]. 

After obtaining a vacuum solution outside the object, we 
can analytically extend it down to all points x (f 7 , in a 
manner described in Refs. El El; because of its l/r n be¬ 
havior, the extended field will not be of uniform accuracy 
very near 7 , and it will diverge precisely on 7 , but it will 
agree with the physical solution everywhere else. 

Since the functions h™„(x',z) depend on z M , we can¬ 
not solve Eq. <[5j) by equating explicit coefficients of e 
without forcing the functional argument z M to be inde¬ 
pendent of e; intuitively, this is the case because the Ein¬ 
stein equation is overdetermined, constrained by the lin¬ 
earized Bianchi identity \7 u (5R l _ liy — \g^g a ^5R a p) = 0, 
which prevents us from freely specifying the center-of- 
mass worldline 7 . So instead of solving Eq. © directly, 
we construct the asymptotic series ([3| by dividing the 
Einstein equation into two parts: field equations that 
(speaking roughly) determine the functionals /i™„(a;;z), 
and constraint equations that determine the functionals’ 
argument z^. 

To accomplish this, we first write Eq. ([ 5 ]) with a gauge 
condition already imposed. Choose a set {E™„} n of 
second-order, linear differential operators that can be ob¬ 
tained from the linearized Ricci tensor via a linear gauge 
transformation; that is, for each tensor there must 
exist a vector field satisfying E™ v [f + C^g] = SR^„[f], 
where £ is a Lie derivative. Define the functionals 
/i” (a;; z) to be solutions to 

E; v [h n ] = S n [h 1 ,...,h n ~ 1 ] (®0 7 ) ( 7 ) 


2 Note that 7 lives in the smooth manifold of the background 
spacetime, not that of the perturbed spacetime, making this 
condition sensible even if the small object is a black hole. See 
Refs. [5l[6llll| for precise formulations in the language of matched 
asymptotic expansions. One works outside the object and im¬ 
poses the matching condition, in lieu of involving the object’s 
interior, to overcome two problems: First, an expansion in the 
limit of small mass and size fails to be sensible at distances ~ e 
from the object (and in the object’s interior); there, the object’s 
own gravity is strong, and the field h^ u it contributes to g^ u 
is comparable to, and has stronger curvature than, the external 
metric g^u- So one should assume an expansion of the form |3j 
only in a region sufficiently far (i.e., at distances e) from the 
object. Second, if we work in the object’s interior, we must spec¬ 
ify its internal structure and composition, and we must design 
different approximations for black holes than for material bodies. 


for arbitrary 7 , subject to the desired boundary and 
matching conditions. For some choices of E there 
may be no solution to Eq. 0 for arbitrary 7 . However, 
if we choose each E™ v to be symmetric hyperbolic, the 
solutions should be guaranteed to exist. So henceforth, 
assume we have made such a choice. 

The sequence of equations 0 form our field equations. 
They are equivalent to the Einstein equation ([5| in a par¬ 
ticular gauge, and hence they take care of a large portion 
of that equation. What is left is a gauge condition, or 
constraint equation, which reads 

Y j e n C; v {h n ) = 0 Or £7), ( 8 ) 

n> 1 

where 

C^[h n ):=SR^[h n ]-E^[h n ]. (9) 

Equations Q and ([ 8 j are the desired division of the 
Einstein equation. By solving Eqs. ([7|, one obtains the 
series 0 as a sum of functionals of two types of free func¬ 
tions: the worldline 7 , which characterizes the object’s 
mean motion, and a set of multipole moments, which 
characterize the object itself and can be defined as ten¬ 
sors on 7 . By enforcing the condition ([ 8 ]), one determines 
the equation of motion Q as well as evolution equations 
for the multipole moments^] These facts about the so¬ 
lutions follow from the algorithmic solution method pre¬ 
sented most thoroughly in Refs. El EEL which I will not 
recapitulate but will draw conclusions from here and be¬ 
low. 

Because each h™ v depends on e (through 7 ), Eq. ^ 
cannot be solved by setting C" [h n ] = 0. Instead, it is 
to be solved by substituting the expansion ([4]) and only 
then solving order by order in e while holding z M (and 
) fixed. This leads to a sequence of equations for 
F“(r; z), in which z M is still held fixed. By holding z M and 

- fixed during this procedure, rather than expanding 
their e dependence, I preserve the particular accelerated 
worldline that satisfies the chosen center-of-mass condi¬ 
tion. Solving the sequence of constraint equations for 
F“(r; z) yields a better and better approximation to the 
equation of motion of that particular worldline, without 
ever expanding the worldline itself. 

In practice, of course, we have access to only a finite 
number of terms Ff in the expansion of the force. If we 
truncate the expansion Q at order n, then we work with 
an approximation zlf to z M ; this approximant satisfies 

^=^e^(r;z n ), (10) 

p—0 


3 For moments of quadrupole order and higher, the Einstein equa¬ 
tion alone does not fully determine the moments’ evolution. The 
evolution of these moments is freely specified, either directly or 
via a specification of the object’s composition. 
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and the nth-order self-consistent approximation to the 
field is a functional of this worldline, given by h liv (x, e) = 
E;Lo eP ^/iv( x ’ z n)- However, to reduce the burden of 
notation, I will typically write expressions in terms of z M 
instead of its approximant z%. 

Thus far I have not referred to the decomposition 
/i M „ = + h^ u with which the paper began. This is a 

decomposition of the extended field, and its singular piece 
h^ v diverges on 7 . A convenient, precise way of choos¬ 
ing the split is described in Refs. HUH El [HJ. With a 
straightforward extension of the analyses in Refs. 13 nu, 
one can show that no matter what the forces F^{t\ z n ) 
in Eq. (10 1 turn out to be, one can always choose the de¬ 
composition of such that the regular field possesses 
the two core properties that (i) Eq. (10) is equivalent to 
the geodesic equation in g^ = g^ v 
a smooth vacuum metric on 7 . 


7c„, and (ii) g^ is 


2. Field and motion through second order 


Let me make this more concrete at first and second 
order. At these orders, the regular fields /i R " satisfy the 
vacuum version of Eq. (|7|): 

J5> M ]=0, (Ha) 

El v [h^} = ^R^\h^\ 1 (lib) 

note that these equations, unlike Eq. |7]), are not re¬ 
stricted to points off 7 . The singular fields satisfy the 
nonvacuum equations 


El[h S1 ] = 8nf^[z\, (12a) 

E%,[h a2 ]=totf%,[z,5m], ( 12 b) 

El v [h Snl } = S 2 R^[h s \h sl ] 

- 25 2 R^h 31 , (i R1 ] (x 7 ), (12c) 


where I have written /i ® 2 as the sum of two parts, h s2 v 


jlV 


the first of which is sourced by a “stress-energy” 
and the latter of which is sourced by nonlinearities. 

The stress-energy source terms I have introduced in 
Eq. (12) are derived quantities defined from the extended 
fields h™ v according to 


T; v ■■= ^E; v [h n ] - s; u [h \..., h"- 1 ]. (13) 

One can loosely, though not precisely [3 HU, think of 
these quantities as the sources of everything in that 
is not generated by nonlinearities. They are supported 
only on 7 , and they are uniquely specified by the mul¬ 
tipole moments that characterize the object. For a non¬ 
spinning object, at the first two orders they (or their 
trace-reversals T^ v := TJf v - \g^ug a0 T^) are given by 


Tp V = / mu^u v 6(x,z)d,T, (14a) 

J -7 

(14b) 


rji2 _ 


8m lxv 8{x , z)dr, 


where S(x, z) := S !yzrg' > is a covariant delta function, m 
is the object’s leading-order mass, 5m^ is a correction 
to the object’s monopole moment, and the bar denotes 
trace-reversal. In addition to the equation of motion for 
7 , the constraint equation ^ dictates the specific forms 
of T" determining a (gauge-dependent) expression for 
and determining that m is constant. Note that T* 
is simply the stress-energy of a point particle of mass 
m moving on 7 ; in this sense, the small object can be 
thought of as a point mass. 

Unlike Eq. |7|, the equations (111 and (12) include 7 
in their domain, with the exception of Eq. (12cI. Equa¬ 
tion ( 12 c) is restricted to points off 7 because the non¬ 
linear source term in it is generically ill defined as a dis¬ 
tribution on any region including 7 . In the same way, 
the total first-order field h 3 = h 3 l + h^l satisfies a field 
equation with a distributional source, 


El„[h 1 ]=8irT^[z], (15) 

but the total second-order physical field hj lL , = h s2 + 
h Snl + /i R2 does not, instead satisfying an equation with 
a pointwise source, 

El v [h 2 ] = -6 2 R^[h\h 1 } (xti). (16) 

However, by moving all curvature terms to the left-hand 
side, we can write the field equation in the distributional 
form 


El„[h 2 ]+6 2 R^[h 1 ,h 1 ] = 87 tT 2 v . (17) 


The left-hand side is well defined as a distribution be¬ 
cause the non-distributional singularities in 5R^ v [h 2 ] 
cancel those in 5 2 R fl ^[h 1 ,h 1 ]. We can also combine the 
first- and second-order equations to fomfl 

di? At!/ [e/i 1 +e 2 ft 2 ] + e 2 5 2 R^h 1 , h 1 ] 

= 87ref^ + 87re 2 f 2 v + 0(e 3 ). (20) 


All of the above equations involve functionals of 7 . The 
final ingredient in the second-order-accurate approxima¬ 
tion is the equation of motion (10 1 . It can be written 
in t erms of the regular field as Eq. (42). As shown in 
Sec. Ill A1 this is equivalent to the geodesic equation © 
in the smooth vacuum metric g^ u . 


4 More generally, we may write the complete Einstein equation as 
the distributional equation 


E e " (^>3 - ( 18 ) 


or simply 


R u 


E 




= 8 ttE< 


( 19 ) 


where all tensors are understood to live on the background. 
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3. Gauge choices and practical implementations 


B. Gralla-Wald 


The structure of the approximation becomes more 
transparent in a particular gauge. In the Lorenz-gauge 
formulation, we have E% v [h n ] = + v p h^ p ) 

and C™ v = for all n, and Eq. Q can be re¬ 
placed with Yh n = 0. Near 7 , the first-order 

singular field (and hence h 1 ) behaves as 

Oivrj 

h-ll = — (g^ v + 2 u^u„)+O{r 0 ), (21) 

where r is a geodesic spatial distance (proper in g M „) from 
7 . The second-order held behaves as 


I next consider the Gralla-Wald approximation, named 
after the authors of Refs. mi- This approximation con¬ 
sists of a strict Taylor expansion, 

VO, e) = z n Kv( x )\ (23) 

n> 0 

here, unlike Eq. ([ 3 ]) , the coefficients do not depend on 
e. Since the metric depends on the object’s motion, this 
expansion of requires an expansion of the worldline 
itself, 

z»{s,e) = Y J t n K{s), (24) 

n>0 


h%= 5 -^+O(r 0 ), ( 22 a) 

C~5r + °( r-1 )> ( 22b ) 


where Sm is given explicitly by Eq. (133) of Ref. 124] . 
Later sections will refer to these local forms to elucidate 
several key ideas. For now, note their essential char¬ 
acteristic: they diverge on the object’s self-accelerated 
center-of-mass worldline 7 , rather than on a background 
geodesic. This distinguishes them from the singular field 
of the Gralla-Wald approximation described below. 


Let me summarize. The self-consistent approximation 
described in this section is based on splitting the Einstein 
equation into a sequence of hyperbolic equations together 
with a constraint. The constraint determines the equa¬ 
tion of motion of the object’s center-of-mass worldline 7 , 
and it constrains the evolution of the object’s multipole 
moments, which are tensors on 7 . From the analysis in 
Sec. II of Ref. [7], one can expect that at a formal level, 
this is the only required input from the constraint: if the 
Cauchy data satisfies the constraint, then the constraint 
should be preserved by the coupled system comprising 
the hyperbolic equations, the equation of motion, and 
the evolution equations for the multipole moments. In 
practice, the set of coupled equations can be solved nu¬ 
merically with a puncture scheme [i] co mum, in which 
one uses residual field variables := h™ — h in 


place of Here is a puncture field that mimics 
/i®" near 7 , such that h and d a h Jf" = d a h^ff on 
7 , allowing one to use h in the equation of motion (42). 


No serious study has been made of constraint viola¬ 
tion or numerical stability in these puncture schemes, 
particularly in the wide class of gauges considered here, 
but such issues are not of essential interest in this paper. 
The question of interest here is “Suppose two researchers 
obtain asymptotic solutions of the form ([ 3 ]), in whatever 
manner, in two different gauges. How are their two so¬ 
lutions related?” This is the question addressed in the 
later sections of this paper. 


where s is a parameter on the worldline. The leading- 
order term, z(f(s) := Zq(s) = z^(s, 0), is the coordinate 
description of a worldline 70 , and the subleading terms 
z £>0 are vectors living on 70 that describe the deviation 
of the perturbed worldline 7 from 70 . Hence, in this 
approximation, one treats the effect of the self-force as 
a small perturbation of the worldline itself, rather than 
as a small perturbation of the equation of motion for . 
Hence, instead of seeking an equation of motion a la Q 
for z ^, one seeks evolution equations for each quantity 


D 2 % 

d T o 


(25) 


where To is proper time (as measured in g^) on 70 . 

It is easy to intuit that one can derive expansions of 
the form (p3| and ([25| from the self-consistent approx¬ 


imation simply by substituting the expansion (24) into 
the field ([ 3 ]) and equation of motion Q . This intuition is 
borne out by the explicit expansion procedure of Sec. m 
On the other hand, one cannot go the other way and ob¬ 
tain the self-consistent from the Gralla-Wald approxima¬ 
tion (although Ref. [5J contained a heuristic argument for 
how to make that leap at first order). 

However, one can certainly work entirely in the Gralla- 
Wald expansion without ever relating it to the self- 
consistent one. The general procedure is very similar 
to that described in Sec. |II A 1[ because of the similarity, 
I shall forgo a general description and simply emphasize 
key points and concrete equations through second order. 
Starting from the expansion (23), one obtains results very 


similar to those of the self-consistent expansion, but with 
a significant modification due to the different treatment 
of the worldline. In the self-consistent expansion, the 
fields h™ v are functionals of 7 and of a set of multipole 
moments on 7 ; in the Gralla-Wald expansion, the fields 
h™ are functionals of 70 , not of 7 , and of a different set 
of multipole moments on 70 . The multipole moments are 
modified by the deviation vectors z!^ >0 , which themselves 
correspond to mass dipole moments and which alter the 
other moments. 
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This is most easily understood by starting at the end. 


Consider Eq. (21), which shows that h^ v behaves like a 
Coulomb field; in flat space it would read oc > 

where ( t,x l ) are inertial Cartesian coordinates. If we 
treat the worldline as in Eq. (24), then . f m it becomes 

- \X Z | 


em _|_ e mzij^x 3 -z J a ) 


l^-zjl 


\ xi ~ z 'o \ 3 


-0(e 3 ). In words, the held in the 


self-consistent expansion, , ,f”T, , diverges at the object’s 

\ x z I 

perturbed center-of-mass position the Held in the 


Gralla-Wald expansion, 


+ 


e z mz\j ( x 3 —Zq) 


+ 0(e 3 ), 


l^-^ol ' \ xi ~ z o \ 3 
diverges at the object’s zeroth-order position Zq, and the 
deviation of the object away from that position alters the 
functional form of the metric, introducing a mass dipole 
mzij. Or as 4-vectors, 


moment Mj 


\ / = mz?,, 


(26) 


where z^ ± = {g^ v + UqUo^z^ is the projection of z\ 


orthogonal to 70 ; here Uq := with r 0 being proper 


dr 0 ' 


time (as measured in g^) on 70 . Equation (26) is in fact 
how we define the object’s deviation from 70 : z^ ± := 
M^/m. 

Concretely, the local behavior of the Helds in the 
Lorenz gauge is no longer given by Eqs. (21) and (22) 
but by 


/cji 2 TO. 

C = —(9, 


h s2 = 


r 

2mzux 


v 2 m 0ai u 0 
( 9n" 


O(r 0 ), 


2u 0 ,u 0 „)+ S ^. 


0(r 


(Snl 




(27a) 

°), 

(27b) 

(27c) 


for the singular field. The stress-energy terms are defined 
in analogy with Eq. (131, and they (in their trace-reversed 


form) are given by distributions on 70 , 


T^= / mu 0lJl u 0l ,5{x,zo)dTo, 


(31a) 


rri2 
^ fllS 


7 7o 


+ muo^uovZ^Vy 5(x, z 0 )dT 0 . 


(31b) 


Just as in the self-consistent approximation, we can 
also write the Einstein equation in distributional form, 


5R fll ,[eh 1 -\-e 2 h 2 ) + e 2 5 2 R^[h 1 , fi 1 ] 


= 8neT^ + 8ne 2 T 2 „ + 0(e 3 ). 


(32) 


The second-order approximation is completed by the 


evolution equations (25) for z%. As in the self-consistent 


case, for a nonspinning object with vanishing quadrupole 
moment, the evolution equations can be written purely 
in terms of the regular field: in place of Eq. (42), we 
here have the evolution equations (541, (60), and (|64|, 


with Eqs. (59b) and (65b). The Riemann terms in these 


equations of motion are geodesic-deviation terms; they 
correspond to the fact that even in the absence of a force, 
two neighbouring curves z M and z q will deviate from one 
another due to the background curvature. 


III. FROM SELF-CONSISTENT TO 
GRALLA-WALD 


where x' 1 are local spatial coordinates centered at 70 and 
r is a geodesic spatial distance from 70 . The monopole 
moment 5m is modified by the presence of zu, and it 
is now given explicitly by Eq. (145) of Ref. [23] 

In any gauge, the field satisfies equations analogous to 
Eqs. ([15]), @, ([IT]), and 


5R^[h 1 }=8nf 3 v [zo], (28a) 

5R^[h 2 ]=-5 2 R^[h\h 1 ] (a^ 7 o) (28b) 


I now show how to obtain the Gralla-Wald approx¬ 
imation from the self-consistent one. This means two 
things: how to obtain equations of motion for z£ given 
an equation of motion for z^, and how to substitute the 
expansion of z M into functionals such as h™ v [z\. 

Two main tools are required in these derivations: the 
ordinary Lie derivative, £, and a variant of it, £, to be 
described below. 


for the full field, 


5R flv [h m ] = 0, 

(29a) 

6R tlv [h R2 ] = -6 2 R liV [h m ,h m ] 

(29b) 

for the regular field, and 


5R^[h sl ] = 8nT 3 „[z 0 ], 

(30a) 

5R^[h s2 ] = 8nT 2 u [zo, 5m, zi], 

(30b) 

5R^[h Snl ) = -5 2 R^[h sl ,h sl ] 


- 2 5 2 R^[h sl , fi R1 ] {x 70 ) 

(30c) 


A. Expanded forms of the equation of motion 


The difference between the self-consistent and Gralla- 
Wald approximations is most conspicuous in their equa¬ 
tions of motion. In this section, I delineate that differ¬ 
ence. I start from the geodesic equation in a smooth 
metric g^ V1 and in Sec. Ill A 1 I derive its self-consistent 
expansion, given by Eq. (42) below. In Sec. Ill A 2 I 


start from the self-consistent result and by expanding 
the worldline, I derive the Gralla-Wald expansion of the 
geodesic e quati on, g iven by Eqs. ( |54| , (60), and (641, 
with Eqs. (59b) and (65b). 




































1. Self-consistent expansion 


Generically, the geodesic equation ([2]) reads 
d ~ 

— + T^zO = kzr (33) 

where s is a potentially non-affine parameter on z^, := 
^- is its tangent vector field, f(j p is the Christoffel sym¬ 
bol corresponding to g and k(s) := ^ In ^/—g ll ^z fJj z 1 '. 
Here for simplicity I have dropped the unknown 0(e 3 ) 
error term in Eq. ([2]); this also serves to emphasize that 
the derivation applies to the more general problem of 
expanding the exact geodesic equation in any smooth 
metric. 

Now, if we write the metric as the sum of two pieces, 
+ h R u , and if we take s = r, the proper 
time on z ^ as measured in g and if we rewrite the 
geodesic equation in terms of covariant derivatives com¬ 
patible with <? p „, we find 


£>V 1 

dr 2 


= -c» vp u v u p + ~ KU>1 = : e )> 


(34) 


where 

C a M := fr^ 7 = \g ab (2/i^. 7) ^ hf r>s ) (35a) 


is the difference between the Christoffel symbol associ¬ 
ated with and the one associated with g Ml/ . With r 
as a parameter, k becomes 


k(r) = 


£ 1-hKuW 


\J 1 - 


(36) 


While z p is a geodesic of Eq. (34) tells us it is accel¬ 
erated in g M „. 

So far no approximation has been made; Eq. (34) is 
exact. If we now expand and k in powers of 

we find 


F* “ ^(n aS h RaS \f n h R — 

- hp^u 01 D ^ 2 u 7 - \ ) h^ 1 . s u a u l3 * u l u & 

- \h%h%. s u a uWvfiu v + 0[(h R ) 3 }. (37) 


Here already is the spirit of the self-consistent expansion: 
by expanding in powers of h R v rather than powers of e, I 
avoid expanding z p {t, e) itself. However, Eq. (37) is com- 
plicated by the fact that the acceleration appears 

in a nontrivial way on the right-hand side. To disen¬ 
tangle it, I now explicitly substitute the self-consistent 


expansion QJ^] After simple rearrangements, one finds 

K = 0, (38) 

F T = -\p aS (2 h^ n) - h%. s ) u?u\ (39) 

n = - l -p^(2h%. n) hf r ,yv? 

+ \p ap hf s (2 h R lp. l} - hf rt5 ) u?u\ (40) 

and hence, 

r yj4- =eF?(T;z)+e 2 F£(T;z) + 0(e 3 ), (41) 


where P ap := g aM + u a u^ projects orthogonally to iF. 
This result may also be written in the more compact form 


2 _a 


D 2 z‘ 
dr 2 


1 


= (<?/ - O (2 h%. n - h R TtS )u^ 


2 

0(e 3 ). 


(42) 


Equation (41) is the self-consistent approximation to 


the geodesic equation. It is defined by expanding Eq. (34) 


in powers of e while holding the solution z^(r,e) to that 
equation fixed. 


2. Gralla-Wald expansion 


In the previous section, I expanded the geodesic equa¬ 
tion while holding ^(r, e) fixed. I now expand £ m (t, e) 
as well. This procedure yields a sequence of equations 
for the vectors z p that measure the deviation away from 
the zeroth-order worldline in the Gralla-Wald approxi¬ 
mation. 

In previous work 03 m nn, this expansion has al¬ 
ways been performed after choosing some set of coor¬ 
dinates and some parameter s (which may differ from r) 
on the worldline. Here I wish to develop a more gen¬ 
eral and more geometrical picture. First, let me set the 
stage. Consider the two smooth metrics and g p „ on 
a smooth manifold Now consider a family of world¬ 
lines 7 e : s £ 1R 1 — > 7 e (s) C A4 with parameter s and 
coordinates z p (s) = € (s)). Each family member sat¬ 

isfies the self-consistent equation of motion (42), which I 
now write as 


< = Ffi(re), (43) 

where r e = r e (s) is proper time on %, a p := is the 
acceleration of and Ffi(r e ) is given by the right-hand 


r 7")2 ~ a 

0 One could instead assume that 2 possesses an expansion in 

powers of 

6 In Sec. |lV| I will put the perturbed and background metric on 
different manifolds, but here it is simplest to consider a single 
manifold. 
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side of Eq. (|42j) . The family generates a two-dimensional 
timelike surface S C Ad with coordinates (s,e), embed¬ 
ded in M according to x^(s, e) = z£(s), which I will also 
write simply as z^(s,e). 

Now if we consider an expansion in powers of e, we 
immediately run up against two apparent ambiguities. 
We can write the expansion as 


^(s,e) = ztf(s) + ez^{s) + e 2 z£(s) + 0(e 3 ), (44) 


where 


1 B n z^ 
Us) = — e " 
lU n! de n 


M). 


The first ambiguity in this expansion stems from the fact 
that we are expanding not a vector but a scalar field, 
equal to the /zth coordinate field on the worldline. Be¬ 
cause of this, the expansion depends on the choice of 
coordinates. How does this play out? The zeroth-order 
term in the expansion is Zg(s ) = z /J, (s, 0 ), which, while it 
is a set of coordinates along a curve 70 , is invariant in the 
sense of uniquely identifying that curve. The first-order 
term, ^f(s), is a derivative along a curve (of increasing 
e and fixed s) in S. This derivative is evaluated on 70 , 
making it a vector on 70 ; again, this is a covariant quan¬ 
tity. But at second order and beyond, the deviations 
z£ lose their vectorial character: unlike the first deriva¬ 
tive along a curve, second and higher derivatives are not 
vectors. The function z fJ, (s , e) describes a curve in a par¬ 
ticular set of coordinates, and the corrections z% depend 
on that choice of coordinates. 

The second ambiguity in the expansion stems from the 
fact that we are expanding at fixed s. Geometrically, we 
are taking the limit e —> 0 along curves of fixed s in 
S. Under a reparametrization s —¥ s'(s, e), the terms 
Zn>o are altered, becoming derivatives along a differ¬ 
ent curve, a curve of increasing e along which s', not 
s, is constant. For a general “small” reparametrization 
s -7 s' = s + O(e), this is a type of gauge freedom, 
similar to but distinct from the usual gauge freedom of 
perturbation theory: rather than a small transformation 
of the extrinsic coordinates x M , it is a small transforma¬ 
tion of the intrinsic coordinates (s, e) on S. The effect of 
reparametrization is displayed explicitly in Appendix [A] 

To avoid a mire of coordinate- and parametrization- 
dependent results, I will isolate the coordinate- and 
parametrization-independent content of each z£. These 
are the quantities I will derive evolution equations for, 
and as will become apparent in Sec. Ill B[ they are the 
only quantities necessary for a practical Gralla-Wald ap¬ 
proximation. 

We can get at these desirable quantities through more 
obviously geometrical ones. First note that the surface S 
can be generated by the two vector fields i M (s, e) = dz>1 


ds 


and zU(s, e) := The field zP describes the deviation 
between neighbouring curves and 7 e +de, and the first- 


order term iff is simply its restriction to 70 , 


z? = v» I 


'7o 


(46) 


The direction of this quantity plainly depends on the 
choice of parameter s. As shown in Appendix [Aj with 
a reparametrization we can freely adjust the piece of 
zf that lies parallel to 70 , but we cannot alter the 
piece perpendicular to 70 . Hence, a covariant and 
parametrization-invariant measure of deviation is 


z i± ■— r 


(47) 


No matter the coordinate system and parametrization in 
which the expansion (44) is performed, the coordinate- 
and parametrization-invariant piece of £f can be picked 


(45) out using z± ± = Pq v z\, where P^ := + ztg 1 


'0 ■ 


For the second-order term, we can construct a “devi¬ 
ation of the deviation” from a derivative of v 1 '. First 
consider the vector 


2 de 

and its restriction to 70 , 


1 Dv a 1 , 0 

w“ := —— 7 — = -v 1 V pv 


^2N w ho- 


(48) 


(49) 


The subscript “N” stands for “normal”: is the 


second-order term in the expansion (44) when that 


expansion is performed in locally inertial (i.e., nor¬ 
mal) coordinates centered on 70 . In those coordinates, 

\v p S7pv a \ l0 = \v?dpv a \ l0 = If^Cs.O) = z%{s). 
Unfortunately, as shown in Appendix TAJ not only z^ 
but also zf NJ _ = Pq depends on one’s choice of 
parametrization. 

Fortunately, we can easily construct a related quantity 
that does not depend on parametrization. This is done 
by projecting out all information about the flow parallel 
to 7 . Replacing zA with its orthogonal projection z;f := 
P ll v v v , and then taking the orthogonal piece of the result, 
we obtain a quantity 


* 2 “ “ \p a y ± Vpvl | 7o (50) 


that Appendix [A| shows is paramatrization-independent. 
z 2 -j- is simply related to according to 


z 2% ~ ^2N_L + -U||^i_l> 


(51) 


where z x \\ := ffzz 0At , zzf := and u% ± := P 0 M 


0 ^“1 


UT = 


Dz 


dr 0 


If we construct z% according to Eq. (44) in some 


coordinate system and with some choice of parameter s, 
then we can extract the covariant and parametrization- 
invariant content of zf using the relation 


° 2 J 


= P, 


0 P 


+ + 


-- i||ZU 


(52) 


v-/ 1 T~\ ^ Q 

which follows from z^ = h dt 2 


ipa 

2 ^1 


and Eq. (51). 


z i± and z^ will be my first- and second-order mea¬ 
sures of the deviation of the accelerated worldline from 
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the zeroth-order geodesic, 
seek evolution equations for. In Sec. |IIIB 


They are the quantities I 
I will show 


that they suffice to obtain the Gralla-Wald approxima¬ 
tion (23) in any coordinate system, despite the coordinate 


dependence of Eq. (44). 


Evolution equations for z x± and can be derived in 
several ways, but here I wish to use a technique that will 
also apply directly to the discussion of gauge transforma¬ 
tions in later sections. And as will become clear in those 
later sections, that means performing an expansion along 
a flow generated by a vector field. Let ip e be the diffeo- 
morphism describing the flow generated by v and let A 
be a smooth tensor field on S (suppressing indices on A 
for compactness). The expansion of A along a flow line 
beginning at Zq ( s ) reads 


(<ptA)(zo(8)) = (e^A)(z 0 (s)) 

= 4l(~o( s )) + t£vA(zo(s)) 
1 


^ClA(z 0 (s)) 


0{e d 


(53a) 


(53b) 


where (p* is the pullback. 


We want to apply this expansion to the equation (431, 


and from sequential orders obtain the evolution equations 
for Zq, z x± , and z!^. For convenience in performing the 
expansion, I choose the parameter s such that it reduces 
to To on 7 o. This does not imply any loss of general¬ 
ity, since it does not restrict the direction of iF; it is 
analogous to choosing a particular set of background co¬ 
ordinates in perturbation theory, which does not restrict 
the choice of gauge. 

Now, the zeroth-order term in the expansion is simply 
a M (s, 0) = F M (s,0). Since F^(s, 0) = 0, this gives us the 
geodesic equation in the background metric, 


D 


2„P 


dT 0 


= 0 . 


The first-order term is 


A-a M (s, 0 ) = £„F(s, 0 ). 


(54) 


(55) 


To evaluate the left-hand side, I write the acceleration 
explicitly as 


a M (s, e) := 


D 2 z» _ /ds \ 2 r; , 
dr 2 \drJ 


= ~ 5 — Z — KZ‘ 




(56) 


where z M := 


dz* 
ds ’ 


Dz* 


and k(s) := 


£ln I then apply Eqs. (|B2]>, (|B3|), p 6 | ), 

and the Ricci identity. The result is 


C. 


,a M =(tV.„-(iVJ 

i 7o \ ' 170 


D 2 %± 

dr ,0 


R» 


ai'fU'U '0 ^1_L^0 ' 


(57a) 

(57b) 


To evaluate the right-hand side of Eq. (551, I first sub¬ 


stitute the expansion (23) into the right-hand side of 


Eq. (|42j) to get the force in the form 

F tl (s,e,'y e ) = e/f(s,e) + e 2 f£(s,e) +0(e 3 ), (58) 


where f x is given by Eq. (39) with the replacement 
h ^l -t and /£ by Eq. (40]) with 


uKn 

_ ’ 

The forces are ordinary tensors on j e , whereas the 

forces Fi‘' are tensor-valued functionals. I leave it to the 


next section to show how the expansion (231 is explic¬ 
itly obtained from the coefficients h^l L (x; z ^) of the self- 
consistent expansion. So now noting that C v e = 1, we 
find 


CvF* I =-\p< 


/ 


'7o 


olRI _ tRl 
za 5 / 3;7 n p ~/-,8 


ft 7 

u o u o 


=: Ff. 


Putting these results together, we get 

— fz a _ p“ „ v v 

— r l -ft iifiv z l±. ^0 • 


D 2 Z i“l 


dr 0 


(59a) 

(59b) 


(60) 


Note that this result is independent of the choice of pa¬ 
rameter s. For any choice, the expansion of the equation 
of motion along the flow of v^ only determines the evo¬ 
lution of z °the piece of z x parallel to itg, which can 
be arbitrarily altered by a reparametrization, is not de¬ 
termined by the equation of motion. 


The second-order term in the equation of motion (43) 


is 


l -C 2 X{s,Q)= l -C 2 v F^ s ,Q). 


(61) 


To evaluate and simplify the left-hand side, I again begin 
with Eq. (56) and then repeatedly apply Eqs. (B2), (B3), 
(B 6 ), and the Ricci identity. The result is 


■zCla a \ 


D 2 %n± 

dr o 


p op rp rfv 


- 2 P 0 “ p R p z\U$ + Mo® 1 ||'^ 1| || 

+ w“ai|| + UqF^ui^ + 2F“u 1 || 


(62) 


where Ui|| == u 0 p.^ and a x \\ := Uq^RAA. 


On the left- 


hand side of Eq. (611 we have 
1 


-£^| 7 o =/ 2 7 S ,0)+A/r( S ,0) 

= /2 0 , 0 ) + ^l/pa( S ,°) 

-w.L- 


(63a) 


(63b) 


Combining Eqs. (|62|) and (|63| in Eq. (|61|), and rewriting 
the 
at 


the result in terms of the variables z xl _ and , we arrive 


D 2 Z% ± 

= F? - PopR p tf. 


v (< z 2t M 0 A 2^1_|7l_L U 6 


) 


+ 2 P^RP^zilu^z 1 ^. 


(64) 
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The force appearing in this equation of motion is 
given by 


F% := /£(«, 0 ) + z? ± fUs, 0 ) - u$F lv i%. 


1_L 


\PS' (: 

W ( 


R2 \ 
v p<r;X ,L (j\;p J 


2 h™. x - h 


*o “o 


R1 


— h 


R1 


pa-XS — 'VA;p« ) “0 U 0~1_L 


er A 


(65a) 


- { 2h fa-x ~ ^U) (5«ij_««o + ^o“ / “i± u o ) ) 
+ (2^ ; a - ^A;p) « (65b) 


Th e seq uence of eq uations (541, (60), and (64), with 
Eqs. (59b) and (65b), are the Gralla-Wald expansion 
of the geodesic equation ((421). They are covariant 


and reparametrization-invariant evolution equations for 
the deviation of the accelerated worldline from the 
geodesic 70 . As such, they represent a refinement of 
the coordinate- and parametrization-specific expansion 
in Ref. [ 8 ] and of the covariant but parametrization- 
dependent one in EH- Other covariant expansion meth¬ 
ods have been used for somewhat similar purposes [29], 
but the analysis here has the advantage of clearly excis¬ 
ing the parametrization dependence of the Gralla-Wald 
expansion. More importantly for the purposes of this pa¬ 
per, its use of Lie derivatives exactly parallels their use in 
gauge transformations, which in later sections will help 
to clarify an essential notion: the gauge-dependence of 
the worldline. 


B. Expansions of the metric and field equations 


The preceding section illuminated the relationship be¬ 
tween the equations of motion in our two formalisms. But 
of course, these equations of motion are largely meaning¬ 
less in themselves; to fully understand the relationship 
between the formalisms, one must know the relationship 
between their metric perturbations h™ v and In this 
section, I show explicitly how the metric perturbations 
(and their held equations) in the Gralla-Wald approxi¬ 
mation can be obtained from the self-consistent approx¬ 
imation by expanding the worldline. 


1. Expansion of functionals of 


written in the forrrQ 


A(x;z) = / B{x, z{s))J-g^v-zt 1 'z v \ds, 
J -y 


( 66 ) 


where B(x,z) is a bitensor that behaves as a scalar at 
z tJ ’(s) and as a tensor of rank (jp,q) at x and primed 
indices indicate evaluation at x M = z M (s). (See Ref. [TO] 
for a pedagogical introduction to bitensors.) Concrete 
examples of such a quantity are T^ v (x] z) and h^^x] z), 
shown below in Eqs. ( |76| ) and (81). We can also write the 
functional as A(x\z) = J B(x, z(s))ds , where B(x, z) := 

B(x, z(s))y/—g ll ' v 'i ll ’z v ' is a tensor at x M and a scalar 
density with respect to reparametrizations of 7 . 

To facilitate the expansion of A, I introduce a Lie 
derivative £ that acts on the functional dependence on 


£fA(x;z) :=—A(x;z + A £)| a=0 (67a) 

= f C^B(x,z{s))ds. (67b) 


This Lie derivative is closely related to the ordinary one 
that acts on the dependence oni' 1 . £ drags fields along 
a flow of the held point x M relative to the worldline 
£ drags helds along a how of the worldline relative to 
the held point x M . 

An expansion of A( x\ z(s, e)) in the limit e —^ 0 is an 
expansion along a how of increasing e in S. We can write 
this in terms of £ as 


A(x; z e ) = ^2 e n S n A(x;z 0 ), ( 68 ) 

n> 0 


where 

S”A(a:;zo) := —£™A(x;z 0 ). (69) 

n\ 

With simple manipulations, one may reexpress this in 
terms of z f and z(( N to hnd 

£A(x;z 0 )= [ % V ll/ B(x,z 0 )ds, (70) 

J 7o 

fc 2 A(x;z 0 )=l [z^ n \/^B(x,z 0 ) 

J 7o 

+ l%'%'VsVv>B(x,zo)]d8, (71) 

where primed indices now refer to the point Zq(s). 

We can also go one step further and express SA(:r; zq) 
and S 2 A(:r; 20 ) i n terms of the reparametrization- 
invariant quantities z^ ± and z ^. After expressing zf and 


We are ultimately interested in expanding the metric 
perturbations and stress-energy tensor. But hrst con¬ 
sider the expansion of a generic tensor-valued functional 
A(x ; z) of arbitrary rank (p, q) (suppressing indices for 
generality). The functionals of interest to us may be 


7 It is not clear whether the perturbations hdfF 2 can be written 
in such a simple integral form, but they can be written in a 
more complicated, implicit integral form [6], which also allows 
an explicit (if more laborious) expansion of the worldline. 












12 


^ 2 n i R terms of these quantities [using Eq. ( |51| ] and writ¬ 
ing B(x,z) = B(x , z(s))y/—gy u 'z^' z v ', one may elim¬ 
inate all parametrization-dependent content by repeat¬ 
edly applying Eq. (Bll and integrating by parts, utiliz¬ 
ing the Ricci identity, and appealing to the equation of 
motion d60l. The result is 


6 A(x; zq) = 

II 

X 

<1 

B(x,z 0 )dT 0 , 

(72) 






8 2 A(x; zq) = 

■/ 

7 ^ X7 

>B(x, z 0 ) + \z^FiyB(x, Zo) 


J 7o 





+ \ 

'~>u! ^i/' 

X~i-L 

^ fi' ^ is' B (x , Zq ) 

dr 0 . (73) 


Here I have expressed the final results in terms of the pa¬ 
rameter To on 70 , but as discussed in the previous section, 
this choice of parameter on 70 represents no loss of gen¬ 
erality. We see explicitly in these results that z^ ± and z!^ 
are the only required measures of first- and second-order 
deviation. 


We can now plainly identify the content of 6 T“^. It 
contains two types of terms: a term and a 6 term. 
The first simply represents the displacement of the point 
particle mass away from 70 . The second is proportional 
to mug a ; under the usual interpretation of compo¬ 
nents of the stress-energy tensor, this represents the par¬ 
ticle’s density of linear momentum in the direction per¬ 
pendicular to 70 , arising from the fact that the particle 
is not just displaced from 70 , but moving away from it. 


3. Expansion of the first-order field 

Next I turn to the expansion of the metric perturba¬ 
tion. In analogy with Eq. ( [75| , we have 

fyL (%; z 0 ) = fi pi/ (x; z 0 ), (78a) 

h%(x; z 0 , h) = h 2 ^(x-, z 0 ) + 8 h*„{x; z 0 , zQ- (78b) 

Likewise, the singular and regular fields are expanded as 


2. Stress-energy tensor 


We are interested in applying the above expansion to 
two quantities: the metric perturbation and the stress- 
energy tensor. Let us first consider the stress-energy. 

According to Eq. ( 681 , the nth-order stress-energy can 
be expanded as 


Og z) = e n rr(x-, zo) + e" +1 STr (x; z 0 , iq) 


+ 0{e n+z ). 


(74) 


and so the first- and second-order stress-energies in the 
Gralla-Wald formalism are 


frfrzo) = T^(x;zo), 
fg v (x;zo,Sm ) z 1 ) =T^ I '(x-,z 0 , 6 m) 

+ $Tf'(x;z 0 ,z 1 ), 


(75a) 
(75b) 

where T and T% v are given by Eq. (14). 

The unknown term at this stage is 8 T| 1 !y (x; zq, z±). 
To calculate it, I first write T^ v (x]z) in the 
reparametrization-invariant form m 

T^(x-z)=m f g°,(x,z)g l 3 p,{x,z)z a ' z 11 ' 

Jj 

S(x, z)ds 


~ ._, (76) 

vXm'XX z v 

where < 7 “, (x, z) is a parallel propagator from the source 
point x M = z fl (s,e) to the field point x M . One may then 
derive 8 T{ iw from Eq. ( |72[ ). Simplifying the result using 
Eq. (Bl) and the distributional identities V ti >S(x, z) = 
— ff/i'V/i<5(x, z) and g“,.p, 6 (x,z) = 0 [ID], one finds 


8 =m g%,g e p, 

J 7o 


2u^“ ufl5(x,z 0 ) 


- zJ ± glV 7 S(x,z 0 ) 


drn. 


(77) 


h% m = h% m (x-, z 0 ), (79a) 

h % R2 = h% R 2 (x; z 0 ) + 8 h s J R 1 (x; z 0 , zQ. (79b) 


As in the case of the stress-energy, the unknown terms 
here are the 6 terms. I express them concretely by as¬ 
suming the existence of a Green’s function for Eq. (151. 
Given a Green’s function satisfying]^] 


E^G p<J y v , = 8 t (g^g^ - 5(x,x ') (80) 


and some boundary conditions, the first-order self- 
consistent field satisfying the same boundary conditions 
is 


h] iV (x-,z) = J G^v/(x,x')Tf " {x'-,z)dV' (81a) 


G 


[IL'H'lS' 


\J-ga'p’Z 01 ' Z& 


==ds. (81b) 


Applying Eq. (72) and simplifying using Eq. (Bl), we get 


l -2-0) ^l) — ^ I ^i_|_ 

+ G^y^-yu^ Uq zjj^dro. (82) 


Our interpretation of the two terms in ET-f^ informs 
our interpretation of this result: the first term in 6 /i* 
is the contribution to the metric from the particle’s mo¬ 
mentum directed away from 70 , while the second term 
is the contribution from the displacement of the particle 


To make the index structure clear, I write Elf v [G\ as 

puvptrrt 
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away from 70 . In the Lorenz gauge, one can use standard 
local expansions of the Green’s function [TO] to estab¬ 
lish this more directly. Although I omit the details here, 
the result of that calculation E3 is that the G^y^-y 
term in Eq. (821 provides the explicit z]' term in the lo¬ 
cal expression (27bI, while the G^y v ' term provides the 
zf-modfication to Sm^ mentioned below that equation. 


4. Field equations 


Finally, we can now immediately obtain the held equa¬ 
tions (28 1 - (30) of the Gralla-Wald approximation by sub¬ 


stituting the expansions of h™ v and T£ v into the held 


equations (15), (16), and ju} @ of the self-consistent 
approximation and regrouping the results according to 
powers of e. 

The analysis in this section has accomplished two 
things. First, it has shown explicitly how to implement 
an expansion of the worldline in the equations of motion 
and held equations, and the relationship between the so¬ 
lutions. Second, it has shown how to do this in a way 
that is covariant and reparametrization invariant. Since 
the quantities 6 T^ v and S h^ v depend only on the per¬ 
pendicular deviation j_, there is no influence from the 
arbitrarily specihed parallel deviation. 

In the explicit expansions I have shown, I have not used 
the second variation (71) of any quantity; they would hrst 


become involved via terms £ 2 T^„ and S 2 h] lll in the third- 


order Gralla-Wald approximation. However, Eq. (73) 


shows that these quantities will again depend only on 
the covariant and reparametrization-invariant content of 
the second deviation. More generally, it should be the 
case that at any order, from the equation of motion for 
z M , one can obtain equations of motion for covariant, 
parametrization-invariant deviation vectors on Zg, and 
one can then write the expansion ( 68 ) in terms of those 
vectors. 

Perhaps most importantly, as I mentioned in the in¬ 
troduction, these tools and results are not confined to 
relating the self-consistent to the Gralla-Wald approxi¬ 
mation. They can be utilized in any scenario in which 
one wishes to expand the self-consistent worldline around 
some neighbouring (possibly non-geodesic) worldline. 


IV. GAUGE TRANSFORMATIONS IN THE 
SELF-CONSISTENT AND GRALLA-WALD 
APPROXIMATIONS 


In the preceding sections, I have elucidated the two 
approximations’ treatments of perturbative motion and 
the relation between those treatments. I now describe 
the effects of smooth gauge transformations within these 
approximations. As we shall see, the differing treatments 
of the worldline entail differing treatments of gauge trans¬ 
formations. 


Section |IV A opens the discussion with a review 
of gauge freedom in perturbation theory, following 
Refs. [5DH52] . The material is standard, but it estab¬ 
lishes my notation and reminds the reader of the basics, 
which are well known at linear order but less familiar at 
nonlinear orders. Sections IIV Bl and II V Cl then focus in 
turn on the self-consistent and Gralla-Wald approxima¬ 
tions. 

In Sec. |IVDl I discuss some general issues related to 
the freedom of choosing what to hold fixed in the limit 
e —> 0 , which is distinct from the usual gauge freedom 
discussed in Secs. IIV Al IIV Cl 


A. Review of gauge freedom in perturbation theory 

1. Active view 

In perturbation theory, we consider a family of metrics 
or simply g in the absence of a chart. (It will 
be convenient in this section to adopt index-free notation 
for tensors.) The family of metrics lives on a family of 
manifolds Ad e , and a given choice of gauge refers to an 
identification map <j>f : Ado —* Ad e . The identification 
map induces a flow through the family down to the base 
manifold Ado where the background metric g ■= go lives. 

Call the generator of this flow X := yyr- We wish to 
approximate a tensor A at a point q = <f>* (p) £ Ad e in the 
limit of small e. In ordinary (i.e., not singular p~8l l33l 134] 1 
perturbation theory, A is well approximated by a Taylor 
expansion around its value at e = 0, given by 

f n 

W?*A)(p) = (e^A)(p) = J2 -0C£A)(p), (83) 

z ' n! 

n> 0 

where p £ Ad 0 - This expansion depends on the choice of 
gauge X, and we define the nth-order perturbation A;} 
to be 


A*(p) := -{£, n x A )(p). (84) 

n\ 


Now say we work in a different gauge. This corre¬ 
sponds to a different choice of identification map <j>Y ■ 

Ado -A M e or flow generator Y := The approxi¬ 

mation of the tensor A in terms of tensors at the point 
p £ Ado is now given by 

(*r*A)(p) = (e^A)(p), (85) 


and the nth-order perturbation is 


Al(p) := 1(£?-A)(p). 


( 86 ) 


The two expansions (83) and |85|) approximate the ten¬ 
sor’s value at two different points in Ad e , q = <j>*{p) and 
q' = <f>\ (p), but in both cases the terms A* and A), 
are evaluated at the same point p in Ado; the situation 
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FIG. 1. Relationship between points in the perturbed and 
background spacetimes in two different gauges. In the two 
gauges X and Y, the point p £ Mo is identified with the 
two different points q and q' in M t , respectively. Likewise, a 
given point q in M t is identified with the two different points 
p and p' in Mo- Very roughly speaking, the transformation 
generators point from q to q ', and their opposites, 
from p to p'. 


is illustrated in Fig. [l] We now ask how the quantities 
and A\ differ when evaluated at this point p. Their 
difference is 

A A n (p) = )(p) - ^(C n x A)(p)- (87) 

n\ n\ 

The first- and second-order terms are easily expressed in 
the usual form 


AAi=£j 1 Ao, (88a) 

XA 2 = L^Aq + — Aq + (88b) 

where £1 := Y — X and £2 := \[X, Y] are the usual gauge 
vectors. Higher-order terms can be easily worked out. 

For compactness, in later sections I will adopt a more 
familiar notation, writing the gauge transformation as 
A n —► A! n = A n + A A n , where the primed tensor refers 
to the Y gauge and the unprimed to the X gauge. 


2. Passive view 

The notion of gauge described in the previous section 
does not utilize coordinates. However, gauge transfor¬ 
mations are sometimes more conveniently thought of as 
small coordinate transformations. 

Let us introduce a coordinate system x M on Mo and 
use the identification maps <f> ' Y and </> Y to carry the coor¬ 
dinates up to M e . In gauge X, keep the name x M for the 
coordinates on M e , such that the point q = </>; Y (p) has 
a coordinate label x^{q) = x M (p). In gauge Y. use the 
name x'^ for the coordinates on M e , such that the point 
q' = (f>Y has the coordinate label x' M (</) = x M (p). The 
two choices of gauge then become expansions at fixed co¬ 
ordinate values x M and x ,p , respectively. So we may do 


away with the differing identification maps, adopt a fixed 
identification (say, X ), and think of the gauge transfor¬ 
mation as the coordinate transformation x M —> x’^. For 
small e, this is a near-identity transformation. To figure 
out its small-e expansion, first consider the differing co¬ 
ordinate values at the points q and q 1 in the chart x /J ". 
The chart consists of four scalar fields, to which we can 
apply the general expansion ([85]) to find 


x^iq') = x*(q) + etf(x(q)) 


£2 ( x (q)) + o£i (x(q))d v ^(x(q)) 


0(e 3 ), 
(89) 


where I have expressed Y derivatives as £ derivatives us¬ 
ing C x x M = 0, and I have expressed the components on 
the right-hand side as functions of x M (g) using the fact 
that x^{p) = x fJ, (q). Now, rather than an active dif- 
feomorphism on M e , let us consider this as the passive 
change in chart x^(q) 1 —>■ x ,At (x(g)). From the definitions 
above, we have x'^{q’) = x^(q). Rewriting Eq. (89) as 
an equation for x IJ, (q ), we then arrive at 


= x^(q') - e£f(x(g')) 


^(*(«'))-^W5 / )MW«')) 


0(e 3 ). 


(90) 


One may now use this transformation to obtain gauge 
transformation rules. For example, the components of g 
transform as 


f) r r CX c)t^ 


(91a) 


= gp V {x', e) + Ct.g^{x', e) + -C\ g^(x', e) 


0(e 3 ), 


(91b) 


where + e 2 ^ + C^e 3 ). We see that in general, 

the tensor transformation is an expansion along the flow 
generated by Recall that we are now using a single 
identification map, say <, and with that identification, 
component values of g are identical to component val¬ 
ues of {<j>f*g). Also note that Eq. (91) applies even if 
g MJ ,(x', e) is not expanded for small e. But if it is so ex¬ 
panded, then the passive transformation law (91) agrees 
with Eq. (88). 


B. Gauge in the self-consistent approximation 


I now move to the self-consistent expansion. This 
is not an ordinary expansion of the sort described by 
Eq. (83). To understand what sort it is, consider its 


form in a chart x /J ". In a gauge in which no logarithms 
of e appear in the metric, it is plain to see from Eq. ©> 
that the self-consistent expansion has the form of a Tay¬ 
lor expansion at not only fixed coordinate values, but 
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FIG. 2. Gralla-Wald and self-consistent limits, shown super¬ 
imposed. A family of spacetimes is stacked vertically, with e 
running upwards. The top surface is the specific spacetime we 
seek to approximate, with a specific value of e. In that space- 
time, the small object is shown in black. In the limit toward 
e —> 0 , this object shrinks toward zero size, and the strong 
curvature due to the object likewise decreases toward zero. 
The Gralla-Wald approximation shrinks the object down to 
a background geodesic 70 . The self-consistent approximation 
shrinks the object down to a worldline 7 that is accelerated 
in the background. 


FIG. 3. Worldlines in different gauges in the self-consistent 
and Gralla-Wald pictures. The family of perturbed worldlines 
{ 7 e } e generate a surface that runs down through the family of 
manifolds Ad e , terminating at the zeroth-order worldline 70 ; 
when this surface is mapped down to Ado using a particular 
identification map, it becomes the surface S of Sec. |III A 2| 
In gauge X, 7 ,= is identified with 7 * in Ado; in gauge Y , with 
7 ^ . The two gauge-related worldlines in Ado are connected 
by a surface generated by the vector field £ M — e 2 ^2 + 

0(e 3 ). The self-consistent approximation works directly with 
7 ; y or 7 ,? . The Gralla-Wald approximation uses expansions 
of those worldlines around 70 . 


also at a fixed coordinate description of the worldline, 
z 11 = 2 ^( 7 ). Hence, if we imagine working with the ex¬ 
pansion ([ 3 ]) prior to imposing the gauge condition, leav¬ 
ing 7 £ freely specifiable, then when we take the limit 
e —> 0 and move down through our family of spacetimes, 
we shrink the object toward zero size and mass on the 
worldline 7 * := (</> £ Y ) - 1 ( 7 ) in Ado- 

After we impose the gauge condition and fix 7 e ’s de¬ 
pendence on e, the situation becomes more complex. In 
this case, the limit e —> 0 takes 7 ,= to 70 , and it may 
seem that we can no longer imagine the “fixed limit. 
However, this is slightly misleading. The essential as¬ 
pect of the approximation is that for a given value of e, 
the object’s worldline, as identified by the location of the 
metric’s ~ 1 /r" behavior, is at the same coordinate posi¬ 
tion in Ado as in Ad e j^] What is happening here is that 
we map from Ad e to Ado holding y e fixed, shrinking the 
object to zero size on 7 ^ and obtaining the perturbations 
/i" {x;z) that live on Ad 0 . By enforcing the gauge con¬ 
dition, we then ensure that those fields on Ado provide 
an asymptotic solution to the Einstein equation. 

This setup is illustrated in Fig. [2] and it remains true 
even in the case that logarithms appear. 

Under a gauge transformation, we not only change the 
coordinates x M —I but also the coordinate descrip¬ 
tion of the worldline, z M —> z ,At , and that change must 


9 Strictly speaking, lies in AT, the manifold of the effective 
spacetime (g M „, AU). This technicality is required because AT 
may not be be diffeomorphic to Ado in a region including 70 . 


be accounted for in the transformation law for the per¬ 
turbations h™ v . In coordinate-independent language, we 
say that in the two gauges we shrink the object toward 
zero size and mass on two different worldlines in Ado: 
le '■= ) _ 1 (7) in gauge X\ yf := (</>* ) _ 1 (y) in gauge 

Y. Figure [3] shows the relationship between these two 
worldlines. 

In the course of this section, I use these ideas to work 
through transformation laws for the worldline, the metric 
perturbations, the singular and regular fields, and finally, 
tensors that live only on the worldline, such as the self¬ 
force. 


1. Transformation of the worldline 


The transformation of the worldline, y £ 


7 , 


Y 


is an 


immediate consequence of the coordinate transformation 


law (901. Under a gauge transformation, the coordinates 
z^(s, ~e) = 2 :^( 7 e (s)) become 


z'»(s,e) = z»(s,e)-eg(z) 


#(*) - n&{z)d v &{z) 


+ 0(e 3 ), (92) 


where functions of z p are evaluated at z tl (s). 

Since it is easy to get lost in the sign conventions, it is 
worth verifying we have the correct signs. Going back to 
Fig. [I] and Eq. (89), we see that the transformation from 
q to q maps to the left in the figure, in the direction of 
Y — X = G. This is the direction of the transformation 
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between the two points in A4 e . Now look instead at how 
the point q is mapped to two different points in A4q: 
p = {4>^)~ 1 {q) and p' = (</>^) _ 1 (g). These two points 
are related by a map to the right, in the direction of 
X — Y = — £ 1 . Comparing the appearance of p and p’ 
in Fig. [I] to that of 7 ^ and 7 * in Fig. [ 3 J we see that we 
have the correct (minus) sign in front of in Eq. (92). 


2. Transformation of the metric perturbations 


Now we may consider the transformation of the metric 
perturbations. As described above, the expansions in the 
two gauges are performed by writing g M1/ in the coordi¬ 
nates x M or x as g M „(x, e;^) or g'^(x', e; z') and then 
expanding for small e at fixed a; M and or x ,At and 
In terms of h^ := g ^ - g^ and h'^ := g( t! , - g the 
two expansions read 

h^(x,e;z) = ^2 e n h™ u (x; z), (93) 

K^'^’Z') = ^2e n h'^(x'-,z'). (94) 


Transforming between the two coordinate systems, we 
find that g^{x, e; z) and g) it/ (a; / , e; z') are related by 

S'^{x',e;z') = J^ga/ 3 (a :{x'),r,z{z')), (95a) 

= gpv(x', e; z(z')) + Ci-g„ v (x', e; z(z')) 

+ X\z{z')) +0{e 3 ). (95b) 


Now expanding g^ u (x', e; z(z')), we get 


h'^{x',e) = e[hl u (x';z)+C il g lxv (x')\ + e 2 [h^ v {x'\ z) 
■^'{2 9nv(x ) + (x j z) 

+ \£\i9nv( x ')] + ^(e 3 )- (96) 


In Eq. (96) I have not yet expanded z^{z') around z^ 
z ,fl . After we perform that expansion, the perturbat 
reads 


h'^{x' ,e) = + +e 2 [hl„(x';z') 

+ A 2 S 'pv{x') + ( C ^ + £^)h^(x'; z') 

+ + ^( e3 )- ( 97 ) 

where £^ 1 hj liy (x l ; z') = Sfe^ a/; z',£) uses the Lie deriva¬ 
tive introduced in Sec. IIII BI 

The two Lie derivatives in ( C^ + £^ 1 )h 1 )lv {x'\ z') have 
an important interpretation. C ^ moves points relative 
to the worldline, and £^ moves the worldline relative to 
the points; hence, their sum has no net effect on the posi¬ 
tion of points relative to the worldline. More concretely, 
the term C^h 1 ^ introduces a mass dipole moment into 
the metric, £$ >1 hj ll/ introduces an equal but opposite mo¬ 
ment, and in the end b! contains no mass dipole moment 
when written as a functional of the transformed worldline 


2 mSjj_ 

\x i -z i \ ' 


z ' M . The fact that the dipole moment vanishes relative 
to z'^ can be seen explicitly by acting with (£j x + £% 1 ) 
on the leading term in h 1 ^, which is 
in 2m the derivative moves x 1 relative to z l , con¬ 
tributing a mass dipole term 2 ™f ira ‘ with moment 

M l = — while in £^ 2m ^ v the derivative moves z l , 
counteracting the effect of . Note that though their 
effects on the mass dipole moment cancel in this way, 
the two derivatives do not precisely cancel one another, 
because has a different tensorial character at x than 
at z. We can infer the net action of the derivatives from 
the result (D2). 

Equation (971 provides the transformation h—> h'^ v 
of the total perturbation. Comparing it to Eq. (94), we 
read off the transformation h™ 
of e n in the two expansions (93) and (94): 


h'™ v of the coefficients 


z') = hl„(x'-, z') + x'), (98a) 

h%(x'\z!) = h\ v {x'\z') + (£ Sl + £ il )h\ v {x'-,z') 

+ £&9pv(x') + \^9^{x'). (98b) 


Note what might be the most essential characteristic of 
these transformation laws: the metric perturbations in 
any given gauge always diverge on the center-of-mass 
worldline in that particular gauge. So ft” [z\ diverges 
on z^, but h'™[z£] diverges on z W 

With little additional effort, we could extend the anal¬ 
ysis in this section to allow functional gauge generators 
££(x\z). Such transformations are necessary (and natu¬ 
ral) because is often required to transform away from 
a specific metric perturbation h^ v (x m ,z), making ££ de¬ 
pend on h^(x', z) and therefore on z However, to avoid 
overburdening the discussion, I leave this extension to the 
reader. 


3. Transformation of the singular and regular fields 

Though we derived the transformation laws for /i" 
with relative ease, we must take care in doing the same 
for the singular and regular fields h®" and . Sup¬ 
pose that in the X gauge the regular field satisfies the 
“nice” properties that g^ is a vacuum metric and z M is a 
geodesic of that metric. In principle, after transforming 
to the Y gauge, we may split /i(”, into any singular and 
regular piece we like. However, clearly we would like to 
do so in a way that preserves the nice properties of the 
split, such that g'nv — 9^ + h'ph is a vacuum metric and 
z'^ is a geodesic of that metric. 

Appropriate transformation laws can be found by 
recalling that for a smooth metric and geodesic, the 
geodesic equation and vacuum Einstein equation are in¬ 
variant under a generic smooth coordinate transforma¬ 
tion. In our present context, the worldline automatically 
transforms according to the standard transformation law. 
It follows that if we let the effective metric also transform 
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according to the ordinary coordinate-transformation law, 
with = ^^g a p(x{x , ),e;z(z > )), then we 

preserve the nice properties of the regular field. Applying 
this rule, we have 

z') = z') + C^g^x'), (99a) 

C 2 ^! 2 ') = h ™( x 'i z ') + (% + £ ii) h ^(x'; z') 

+ £t 2 9nA x ') + \ E \ 1 9^(x')- (99b) 


At the same time, the perturbations h'™ = h'^f 1 + h'jfff 
must satisfy Eqs. (98). This forces the singular field to 
transform as 


h'^{x'\ z') = h^lix'; z’), (100a) 

ti^(x'iz') = hf v (x'\z') + (£ ?1 + £ il )hll(x'-,z'). 

(100b) 


To reiterate, with these rules, the effective metric in 
all smoothly related gauges is always a smooth solution 
to the vacuum Einstein equation, and the center-of-mass 
worldline is always a geodesic of that effective metric. 
Other transformation laws could also do the trick, but 
these are the most natural. 


4- Governing equations in different gauges 


By design, the transformation laws (92), (99), and 


(1001 ensure that the governing equations of the self- 


consistent approximation are invariant under a smooth 
gauge transformation. 

By this I mean, first and foremost, that in all smoothly 
related gauges, the regular field satisfies the vacuum 
equation 


SR^[eh R1 + e 2 h R2 } + e 2 S 2 R^[h R \ h R1 ] = 0{e 3 ) (101) 


and the worldline satisfies the geodesic equation (42). 


These facts are obvious from the transformation laws for 

(.Rn 

Even beyond the regular field, we can establish that 
the equations of the self-consistent approximation have a 
certain invariant form. For example, we can work out the 
distributional equation satisfied by the full field. Starting 
from Eq. (20), using Eq. (98) to e x pres s h™ u [z\ in terms 
of h'™ v [z'], appealing to Eqs. (C5HC7), and expanding 
Tff v [z] around T™ v [z'], we arrive at 


SR^[eh n + e 2 h' 2 ] + e 2 6 2 R^[ti\ti 1 ] 

= 8nef 3 v [z , ]+8ne 2 f 2 u [z'] 

+ 87re 2 (£ 4l + £^)f^[z') + 0(e 3 ). (102) 

This is simply Eq. ( |20[ ) with a transformed distributional 
source, 


T£,W] = TlM ( 103a ) 

T'llA = T%\A + e 2 0C 6 + (103b) 


Therefore the field equation (20) is satisfied in all gauges, 


although the source must be appropriately transformed 
between them. This change in the stress-energy corre¬ 
sponds precisely to the transformation of the singular 
field in Eq. ( |100[ ). It is examined in detail in Appendix [ d| 
So far we have only seen that that the total field equa¬ 
tions R^[g + h R ] = 0 and R^\g + h] = Y, n hold 

true in any gauge (minding that must be transformed 
along with h™ u ). However, we may like a stronger result: 
that the individual fields: h' R J l \z'] and h!ff v \z’\ satisfy the 
separate field equations (11), (15), and (17). This is non¬ 


trivial because we are not working with ordinary Taylor 
series. In Sec. IIA 1 I defined the functionals h as the 
solutions to Eqs. (|7|), enforcing a particular division of 


hfj ,In writing down Eq. (98) for h'ff [z '\, I have adopted 
another particular division of h! jJV in Eq. (97). A pri¬ 
ori, these two divisions could disagree; eh'^z'} could dif¬ 
fer by o(e) from a solution to an equation of the form 
E'^leh' 1 ] = 8 TreT^[z'], for example. 

To see that Eqs. ©, ©, and ([17]) do hold in all 
gauges, return to the idea that the se f-consistent ap¬ 
proximation corresponds to a Taylor series at fixed z M . 
In accord with that idea, write the operators E™ v in any 
gauge as 


E lA hl [A] = [9 + ^[z\ + 0{e 2 )]\ e=i 


£=0’ 


- 5 2 Rf lv [h 1 [z},h 1 [z\\ 


e=0 


(104a) 


E lA h 2 i z l\ = |^^^[3 + e* 1 [~]+e 2 ^ 2 W + C , (e 3 )] 


(104b) 


where before evaluating the derivatives, the expansion Q 
of the equation of motion is substituted, and the evalua¬ 
tion at e = 0 is performed while holding z^ and fixed. 

With these definitions, let us first establish the desired 
result for h' 3 v [z']\ 


E^[h /1 [z']\=8nf 3 I/ [z'}. 

(105) 

The left-hand side may be rewritten as 


O'V]] = 

(106a) 

= jJRA^iz'}] | 6=0 

(106b) 

II 

(106c) 


where the secon d line is a consequence of SR^lC^g] = 0. 
Equation (105) then follows from E^lh 1 ] = 8nT^ v . 


The desired result for h' 2 v [z'] 


E' 2 [h' 2 [z'}} + 6 2 R li „[h n [z'},h n [z']} = 8 t rf 2 [z'}, (107) 


follows from Eq. (104b) in the same manner [repeating 
the operations that led to Eq. (102)]. And we may like¬ 
wise obtain the desired equations for h'^fflz 1 ]. 
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Let me summarize: If the approximation scheme de¬ 
scribed in Sec. IIA l| is applied in two slightly different co¬ 
ordinate systems related by an equation of the form (901, 
then the terms in the two resulting approximations are re¬ 
lated according to Eqs. (f98l) , (99), and (1001. In all cases, 


the center-of-mass worldline is a geodesic of the effective 
metric, and the effective metric is a vacuum metric. 


5. Transformation of the self-force (and other fields on the 
worldline) 


In the above sections, I described the effect of a gauge 
transformation on quantities evaluated at a given point 
in Ado- However, in self-force theory, we are often inter¬ 
ested in something quite different: the effect of a gauge 
transformation on quantities evaluated specifically on the 
worldline. This differs from the previous case because the 
position of the worldline itself changes under the trans¬ 
formation. 

The correct treatment of this situation can be deduced 
from Fig. [3] Suppose that in gauge A', we calculate some 
tensor A = A 0 + eAi + ... on the worldline 7,5, and in 
gauge Y we calculate the same tensor as A' = A 0 + eA[ + 
... on the worldline 7* . How are the results related? The 
worldline 7^ is mapped to 7* via the transformation 
tp e := (</>5 ) _1 o (ff. We can find the tangent field to this 

map, call it £ **, by writing e C Y = ip* = (</>*)*(<pY ) _1 * = 
e e£x e -e£y_ Expanding the exponentials, we find 

If** = -etf -e 2 ££ + 0(e 3 ). (108) 

We can then expand tensors at 7* around their values at 
7,5 using ip* = e C Y. Hence, the difference between the 
expansions in the two gauges on the two worldlines is 

AA = ip* A 1 - A. (109) 


Here we find the change in the result not only due to the 
change of gauge, A —> A but also due to the change 
in evaluation point, p £ 7,5 = ip e {p) G lY■ The 


comparison is made at the original point on 7,5, and 
that is where the tensor AA is defined. I use a slanted 
A to distinguish this transformation from the ordinary 
transformation A —> A'. 

Equation (1091 becomes clearer when applied to a par¬ 
ticular case. To demonstrate its use, I now work out the 
gauge transformation of the equation of motion a^ = F M . 

On the left-hand side, we have 


Aa» = -eC^aP - e 2 (£ ? X - | J C| i o'*) + 0(e 3 ). (110) 
On the right-hand side, we have 


AF» = eAFf + e 2 [A F% - £ €l (Ff + AFf)] + 0(e 3 ). 

(Ill) 


To evaluate these expressions, I return to the type of 
calculations performed in Sec. Ill A 2 In that section, 


I took Lie derivatives of the self-accelerated equation of 
motion to expand it around its value on a background 
geodesic. Now I expand its value on one self-accelerated 
worldline around its value on another self-accelerated 
worldline, with the two worldlines related by a gauge 
transformation. The role of the surface S in Sec. IIII A 21 
is now played by the shaded surface connecting 7,5 to 
7<T in Fi g- 0 and the role of the vector held v M is now 
played by £ However, the situations are not quite 
identical, partly because the acceleration does not van¬ 
ish on either worldline, but more importantly because 
we cannot easily choose an arbitrary parameter on y5- 
If we parametrize 7,5 with proper time, call it tx , then 
the map ip e (le( T x)) induces a natural parameter s(tx) 
on 7 ^ , defin ed via z '^(s(tx)) = x ti (ip fL ('y e (Tx)))- 

As in Sec. Ill A 2 I rewrite the acceleration on 7^ in 
terms of s. 


ds V 


a>1 = [ dY ^ ~ ****) > 


( 112 ) 


where r is proper time on 7* , and s now refers specifi¬ 
cally to the parameter induced by ip e . With this specific 
parameter, just as we had = CyZ^ = 0, we now 

have 


C t Z ' 


** = = 0, 


and so 






(113) 


(114) 


If we did not use the naturally induced parameter, these 
useful re|ations would not hold. 

Since £ M plays precisely the same role as v M did pre¬ 
viously, we c an a p ply all the results of Sec. [B] Con¬ 
cretely, Eqs. (B2) (B6) all hold true with the replace¬ 


ment v 11 —» £ 

With those identities in hand, along with the expres¬ 
sions for Awe can go ahead and calculate AaP and 
AF M according to Eqs. (110) and (111). Since aT = F^, 
we necessarily get Aa M = AF M , and the explicit calcula¬ 
tions bear this out. Either way, one finds the following 
transformation law for the self-force: 

AF» = -e (P» v g + 

- e 2 [p^c+Wu<u+& tV (lr+ R U < 1U ) 

- 2R IJ 'uZ 1 uiiUr + + u^R u ^ iU ^ 

+ 2-^£i“£i F — 2 ^^ 

- vp (Ira,+4ir u v & Up ) - 2irir«, + 

(115) 


+i? (u^r,, - £?,„)] +0(e 3 ), 


where := £2 — §£i£i-i/> overdots denote covariant dif¬ 
ferentiation with respect to r, and I have used the no¬ 
tation ■= R*aP'r,v£.i uP Qvy and R ^ uilU | ?1 := 
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aM\v ua ii ul fX (for example). The first-order term 
is precisely the standard transformation law derived by 
Barack and Ori m- The second-order term appears here 
for the first time. 


We can also apply the transformation law (109) to any 


other quantity of interest on the worldline. For example, 
applying it to the regular field evaluated on the worldline, 
we get 

Ahp,, = eC^ 1 g )J/L , + e {C^ 2 g^ v — P 0(e )• 

(116) 


Here the + £^)h^p terms in Eq. (99bI are removed 
by the action of ip* , leaving only a transformation of the 
background metric. 

As a final observation in this section, I note that one 
could choose a gauge in which the “accelerated” worldline 
in the background is in fact the geodesic 70. This is 
easily seen from Fig. [3] in place of the vector fields A' or 
Y, simply choose a vector field Z that is tangent to the 
shaded surface connecting j 0 to 7 e . This vector field will 
define an identification map cpf satisfying </>f (7o(s)) = 
7 e (s). Hence, from this perspective the deviation of the 
self-accelerated worldline from the background geodesic is 
pure gauge. This will be illustrated more explicitly in the 
next section. However, I will argue in Sec. [V] that such a 
gauge choice is unacceptable on the physically interesting 
domains of size ~ 1/e; on such a domain, it forces the 
metric perturbation to become large. 


C. Gauge in the Gralla-Wald approximation 

I now consider the effect of a gauge transformation in 
the Gralla-Wald approximation. My treatment mirrors 
that of the previous section. However, since the coeffi¬ 
cients hJpv in the Gralla-Wald expansion are independent 
of e, gauge transformations require somewhat less care. 

Unlike in the self-consistent case, where we imagine 
the object shrinking toward zero size on an accelerated 
worldline in Ado, here the object shrinks toward zero size 
on a background geodesic 70, as illustrated in Fig. [2] 

1. Transformation of the worldline 


where functions of z q are evaluated at Zq(s). Note that 
the zeroth-order worldline is unchanged; the effect of the 
transformation is to alter the deviations relative to that 
worldline. 

These results depend on both the background coordi¬ 
nates and the parametrization of the accelerated world¬ 
line z fJ '(s, e). We are perhaps more interested in the 


1_L 


covariant and parametrization-invariant quantities z‘ 
and ^2t- Projecting out the parallel components in 


Eq. (117b), we get 


4±( s ) = ^i_l( s ) — £i_l(“o)> 


(118) 


where := /. Evaluating Eq. ( |117c| ) in a normal 

coordinate system centered on 70, we get 


4n( s ) = ~2 n( s ) “ £2 (zo) + 2^1 (^o)V^(^o) 




(119) 


This result is covariant but parametrization-dependent. 
We can next apply Eq. (51) to get zf^ = ^n_l + 


»/U 

HI- 


When combined with Eqs. (119) and (117b), this gives 
us (after some simplification) 


_ zu 
" 2 % ~ ~ 2 1 


~2J 


- 

?2J 

-£ 


— c 
- r 0iA 


2_L 


Vatf 


(120a) 

(120b) 


where := + £i|| and ^ := % ± - 

Note that the second formula for the transformation, 
Eq. (120b), involves derivatives of Mq perpendicular to 


70, requiring some extension of Uq off 70; this extension 
is arbitrary, and it can be chosen to make the Lie deriva¬ 
tive term vanish, for example. 


2. Transformation of the metric perturbations 


Now turn to the metric perturbations, 
transformation laws (88), we find 


Applying the 


( x ’ z o) = % 9pu, (121a) 

A ^'pi/(‘ r i“o) = PiiQnv A 2A h ^(x^zq) . 

(121b) 


In the self-consistent case, we had a transformation of 
the entire worldline In the Gralla-Wald case, we in¬ 
stead have transformations of each term zf L in the expan¬ 
sion of zffi They can be found immediately by expanding 
Eq. (92) in powers of e, yielding 


4 M ( S ) = z o(s), (H7a) 

zf(s) = zf(s) - £f(z 0 ), (H7b) 

4 M (s) = ^2 (s) - £2 ( z o) + \fX{z 0 )d^{z 0 ) 

- zf(s)d„tf(z 0 ), (117c) 


All the perturbations h'ff in the new gauge, like all the 
perturbations /(” in the old, diverge on the zeroth-order 
worldline Zq . This is an essential point of difference from 
the self-consistent approximation: In the self-consistent 
approximation, a gauge transformation shifts the curve 
on which the singular field diverges. In the Gralla-Wald 
approximation, the zeroth-order worldline, on which the 
singular field diverges, is invariant. 

Instead of altering the curve on which the fields di¬ 
verge, the gauge transformation alters the singularity on 
that curve, by altering the functions z!f >0 that appear in 
hffff 1 . It effects this change in singularity structure by 
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altering the mass dipole moment. Recall the discussion 
in Sec. |IVB2 


from which we can infer that the term 
£ fl h},„ generates a mass dipole moment — m£\ relative 


-"Cl 

to the worldline on which h 1 ^ diverges, which is z q . In 
the self-consistent case, we began with no mass dipole 
moment in the unprimed gauge, and the dipole moment 
induced by was exactly cancelled by the action of 
£$!, leaving us again with no dipole moment. In the 
Gralla-Wald case, we instead generically begin with a 
mass dipole moment M^, and the gauge transformation 
alters it by an amount = —m^ ± . 

I remind the reader that the first-order deviation in the 
Gralla-Wald approximation is defined as z ^_ L := M^/m. 
Hence, the transformation law AM M = —m^ ± , obtained 
from the transformation of the metric perturbation, pre¬ 
cisely reproduces the transformation law (1181. In fact, 
at any order, the transformation laws for z£ >0 , as given 
at first and second order in Eq. (117), can always be de¬ 


rived directly from the transformation laws for the met¬ 
ric perturbations. In principle, one need never appeal to 
Eq. (92). This becomes important in the case of transfor¬ 


mations that are singular on the worldline, which spoils 
the application of Eq . (|92]) but does not prevent the ap¬ 
plication of Eq. (121) [51 1191 1201 . 


3. Transformation of the singular and regular fields 


We must now apportion Eq. (121) into transformation 


laws for the singular and regular fields. In exact analogy 
with the self-consistent case, we make the natural choice 
of transforming the regular field like an ordinary smooth 
perturbation, which now means 




(122a) 


1 


A r ) £'£ 1 9n 1 ' "b 

(122b) 

This leaves the singular field to transform as 

Ah%{x\ z 0 ) = 0, (123a) 

A h s l(x;z 0 ) =£ i : 1 h s l(x;Zo). (123b) 


As in the self-consistent case, these laws preserve all 
the nice properties of the singular-regular split. 


find that the field equations (28) for /i" are invariant, 


but the stress-energy is altered between gauges. Specifi¬ 
cally, Tfi, = T^, and T^ v = T* v + A T* u , where 


= A 


T 1 . 

1 flV 


(124) 


An explicit expression for TL, can be derived from 


one for C^ 1 T ^ 1 ', given in Eq. (Dl). From that equation, 


we see that the Lie derivative both alters 5m and more 
notably, shifts the mass dipole moment by an amount 
Aj_, in agreement with the discussion above. 


5. Transformation of the self-force (and other fields on 70 ) 

Working out gauge transformations of fields on the 
worldline is simultaneously simpler and subtler than it 
was in the self-consistent case. I first present the simpler 
viewpoint; at the end of the section I turn to the subtle 
point. 

From the simple viewpoint, transformations of tensors 
on 70 are straightforward for the obvious reason that 70 
is unmoved by the transformation. Hence, we can use 
the same transformation laws for tensors at a point on 
70 as at any other point in A4o- 

Let us utilize this perspective in working out the trans¬ 
formation of the self-forces in Eqs. ( |59b| ) and (65b). Since 
the equations of motion (60) and (|64l are valid in any 


gauge, we can easily find the transformations of F^ by 
writing those equations in terms of quantities z[ ± = 
z 1± + Az 1± , z' 2f = z 2 j + A z 2 %, an d F£ = F£ + A F£. 


Substituting Eqs. (118) and (120) and then solving for 
A FE 


yields 


AF 1 = ~ 0, ( 125 ) 


.TCP 


and 




— + R a ^u^AzLuq 


- 2 Pr 


0p 2 


+2 pz p r» ^ 

, (126) 


f. Governing equations in alternative gauges 

By construction, all the governing equations are invari¬ 
ant: the equations of motion (60) and (64) apply in all 
gauges, and in all gauges the regular fields satisfy the 


vacuum Einstein equations (29). 

As in the self-consistent case, the only thing to be 


mindful of is the transformation (123b) of the singular 


field, or put another way, the transformation of the stress- 
energy. Repeating the analysis that led to Eq. (102), we 


where Az^ is given by Eq. (120). 

One could equally well derive these results by dir ectly 
substituting h'^J 1 = 


A into Eqs. (|59b|) and 


(65b) and then using the expressions ( [122^ for A h*£. 

Now to the subtle issue, which is clearly seen from a 
geometrical standpoint. Consider a tensor T defined only 
on the family of worldlines 7 e . Figure [3] shows that the 
relationship between 7 e C A4 e and 70 C A4o is unaltered 
by a gauge transformation: the shaded surface connect¬ 
ing them, call it S *, is gauge independent. It follows 
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that the expansion of a tensor T that lives only on the 
worldline is gauge invariant. Its expansion always reads 


e eCz T I . 

■To. 

To write the quantities in this expansion in any par- 
ticular gauge X, we can relate the situation back to 
the expansion of the worldline in Sec. |III A 2| The sur¬ 
face S described there is a gauge-dependent quantity, 
equal in the gauge X to a surface con¬ 

necting 70 to 7* y . It is genera ted b y the vector field 
v , which in analogy with Eq. (108) is now given by 
v = e(Z — X) + |e 2 [A', Z\ +0(e 3 ). So in gauge X the ex¬ 
pansion reads e tCz T = T 0 + e(Ti + C V T 0 ) + 0(e 2 ), where 
Xi = C\T. In a gauge Y the expansion would read 
e eCz T = t 0 + e(T[ + £ V 'T 0 ) + C>(e 2 ), where T[ = C y T 
and v' = e(Z—Y) + ^e 2 [Y, Z] + 0(e 3 ). Clearly, the results 
in the two gauges must agree, since they both began from 
the gauge-invariant field Z. Put another way, Z provides 
a preferred reference gauge in this scenario, and when 
we examine the expansion of quantities defined only on 
the worldline, we must always transform to this reference 
gauge, with the vector v acting as the gauge vector. 

However, in most circumstances, this subtle point need 
not be minded. In a typical situation, we are faced with 
some quantity written in terms of tensors on 70, and we 
are only interested in whether the value of that quan¬ 
tity is altered when we calculate the tensors in a differ¬ 
ent gauge. We can then go ahead and adopt the simple 
viewpoint described above. 


D. Gauge freedom of another kind 

In general, three kinds of “gauge freedom” exist in 
perturbation theory. In this paper, I have focused ex¬ 
clusively on one of them: the identification of points in 
the perturbed spacetime with points in the background 
spacetime. This is the gauge freedom usually discussed 
in the context of perturbation theory [30U32] ; it is equiv¬ 
alent to expanding out the effect of a small coordinate 
transformation, such that the background is left un¬ 
changed and order-e transformations of background fields 
are treated as alterations of the perturbations to those 
fields. 

A second type of gauge freedom is the freedom to per¬ 
form a coordinate transformation (or diffeomorphism) on 
the background, inducing, via one’s identification map, a 
coordinate transformation in the full spacetime. 

However, prior to these two types of freedoms, there 
is a third type: the choice of the spacetime family itself. 
Given a particular metric g involving a numerically small 
parameter e, one can embed it into different families g € . 
This corresponds to choosing “what to hold fixed” while 
taking the limit e —> 0. As a trivial example, consider 
adding a linear, global mass perturbation h 1 = §^8M 
to a black hole metric g. Such a perturbation can be 
incorporated into the background g via a redefinition of 
the background mass, M —7 M+8M , thereby choosing to 
hold M + 8M , rather than M, fixed while e varies. Since 


the mass of a spacetime cannot be altered by a coordinate 
transformation, this trivial example illustrates that this 
freedom is independent of the usual “small coordinate 
transformations”. 

In the context of self-force theory, this third freedom 
has important consequences. Consider the Gralla-Wald 
approximation. Choosing a family of spacetimes requires 
choosing the family of worldlines z^. As an example, 
suppose we are given a spacetime {g fiv (x), A4) contain¬ 
ing a small object of mass m moving about a nonspinning 
black hole of mass M in a quasicircular orbit. On a short 
timescale the orbit has an approximately constant orbital 
frequency 12* measured by an inertial observer at infin¬ 
ity. We can imagine embedding this spacetime into many 
different families. We might choose a family g p „(x, e; 12), 
satisfying g M „(£, e*; 12) = g^(x) for some e* > 0, in 
which all members share the same orbital frequency 
12 = 12(e*) := 12*. But we might instead choose a family 
g Ml/ (:r,e; 12(e)), satisfying g^(x, e*; 12(e*)) = g^(x), in 
which the orbital frequency varies between family mem¬ 
bers, such that 12(e*) = IT but 12(e 7^ e*) 7^ 12*. 

With all else the same, in both cases we write g^ u {x) 
as an expansion around the same background metric 
g^{x) = g M „(x,0;fl) = g M ^(a:, 0; £1(0)), which in the 
present example is a Schwarzschild metric of mass M. 
But the background geodesic around which we expand 
the worldline differs between the two families. In the 
first case, we expand z M around a circular background 
geodesic of the same frequency, 12*, while in the second 
case we expand it around a circular geodesic of a differ¬ 
ent frequency, 12(0) 7^ 12*. Either way, we describe the 
same physical spacetime (g /il/ (a;), AI), but we do so by 
writing it as an expansion along flows in two different 
families of spacetimes. Because the orbital frequency (as 
measured by an inertial observer at infinity) is a gauge- 
invariant quantity, the two families are mathematically 
distinguishable; the difference between them is not a dif¬ 
ference that can be removed by a small coordinate trans¬ 
formation. 

This example is relevant for many practical calcula¬ 
tions in the literature. For a canonical instance, re¬ 
fer to Refs. [3S (3B]. There the authors first expanded 
the circular orbit about a circular geodesic of the same 
Schwarzschild-coordinate radius r$. This specifies a per¬ 
fectly acceptable family of spacetimes, the members of 
which contain orbits of differing frequencies 12(e). How¬ 
ever, the authors then re-expanded their expressions such 
that the end result is equivalent to expanding the circu¬ 
lar orbit about a circular geodesic of the same orbital 
frequency 12 = 12*, now placing the metric in a new fam¬ 
ily of spacetimes, each member of which contains an or¬ 
bit of frequency 12* H This second choice is more useful 
in most cases, because it enables comparison to other 
formalisms, such as post-Newtonian theory, which can 


10 Of course, one can instead expand at fixed frequency directly, 
without first expanding at fixed orbital radius. 






22 


obtain results at fixed frequencies but, because they do 
not expand around a Schwarzschild background, cannot 
obtain results at fixed Schwarzschild-coordinate radius. 
However, from the perspective of perturbation theory on 
a given background, the “hxed-r 0 ” and “hxed-fl” expan¬ 
sions are equally admissible, and they equally well ap¬ 
proximate the particular spacetime A4). 

For more details of how this type of freedom can play 
out, I refer the reader to Ref. m- 


V. SUMMARY AND DISCUSSION OF 
PRACTICAL ISSUES 


In this paper, my central concern has been elucidating 
the mathematical structure and gauge freedom of approx¬ 
imation schemes involving perturbed motion in general 
relativity. The presentation has been, for the most part, 
quite abstract. However, it has also led to several con¬ 
crete, applicable results, including 


1. a covariant and reparametrization-invariant expan¬ 
sion of an accelerated worldline around a neigh¬ 
bouring geodesic. This led to the second-order 
equation of motion (64) and the worldline-deviation 
terms in the stress-energy tensor <[77j) in the Gralla- 
Wald approximation. Moreover, the methods can 
be applied in any problem in which one is inter¬ 
ested in expanding one worldline about another. 
Such expansions appear in osculating-geodesic and 
two-timescale expansions, amongst others. 

2. gauge transformation laws for the metric pertur¬ 
bation of the self-consistent expansion, given by 
Eq. (98), as well as a general transformation law 


for tensors that live only on the worldline, given by 
Eq. (fl09). 


3. natural transformation rules for the singular 
and r egula r field s, given by Eqs. (99) (100) [or 
Eqs. (122) (123) in the Gralla-Wald approxima¬ 


tion], which preserve the essential properties of 
those fields and leave invariant all the governing 
equations of self-force theory. 

4. a transformation law for the second-order self-force, 


given by Eq. (115) [Eq. (126) in the Gralla-Wald 
approximation]] 

In the remainder of the discussion, I describe some 
practical issues related to these results. 


A. Working in gauges other than Lorenz 

To clarify point 3, I note that to make use of the fact 
that the governing equations are the same in all smoothly 
related gauges, one must first find an effective metric 
satisfying them in a particular gauge. And given that 
multiple effective metrics can satisfy the same governing 


equations, one must realize that under a gauge transfor¬ 
mation, one is referring to the transformation of one’s 
particular choice of effective metric. 

Currently, an effective metric satisfying the desired 
properties through second order has been found only in 
the Lorenz gauge [U [3 |24j. Furthermore, although a lo¬ 
cal expansion of the second-order field has been derived 
in an alternative gauge [5], it is available in a practical 
form only in the Lorenz gauge [24]. 

However, with the results of this paper, in principle one 
can calculate the held in any gau ge sm oothly related to 
Lorenz. One need only apply Eq. (100) to find the correct 
singular held in the non-Lorenz gauge. Of course, this 
requires knowing the gauge generator that brings the 
hrst-order held from Lorenz to non-Lorenz; depending on 
how h is found, obtaining this gauge generator can be 
quite difficult. 

A larger issue is that the restriction to smooth trans¬ 
formations is quite severe. The gauges that are most 
convenient for explicit computations in a black hole back¬ 
ground, the Regge-Wheeler and radiation gauges, are not 
related to Lorenz by smooth transformations. One might 
be able to utilize the simplicity of these gauges by work¬ 
ing in a mixed gauge: adapting the ideas of Ref. 0, one 
might implement a puncture scheme in which the punc¬ 
ture is constructed from the singular held in the Lorenz 
gauge, but the residual held (the numerical variable in the 
scheme) is calculated in any desired gauge. The trans¬ 
formation to such a mixed gauge is not perfectly smooth 
on the worldline, but it is of nonnegative differentiability, 
since it only alters the gauge of the residual held, not the 
puncture. 

If one wishes to work in gauges further removed from 
Lorenz, substantially more analysis is required. In partic¬ 
ular, if the gauge generator is not C l at the worldline, 
or if £/( is n °t C° there, then the transformation law for 
the worldline is not given by Eq. (92). This spoils the 
reasoning used to choose the relation (99) between the 
regular helds in the two gauges. In that case, one must 
derive an alternative to Eq. (92) from hrst principles, and 
from there write new transformation laws for the singular 
and regular helds. This analysis has been carried out at 
hrst order in e in recent formulations of self-force theory 
in singular gauges HMD, but it remains to be done at 
second order. 


B. Gauge and motion on large scales 

Forget, for a moment, the question of which gauge we 
are using, and simply suppose we have computed the 
motion and metric in some gauge. We must next consider 
how the results depend on that choice of gauge. Because 
the worldline is a gauge-dependent quantity, the answer 
is not obvious, and we must do some work to disentangle 
the motion from the gauge. 

First consider the Gralla-Wald case. Here it may seem 
that motion and gauge are entirely indistinguishable, 
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since we can freely adjust each deviation vector z£ with 
a gauge transformation. Nevertheless, we can distinguish 
the two notions. Begin in a gauge where the first-order 
held hr {x\ z o) does not grow large with time; any of the 
typical gauge choices, such as Lorenz or Regger-Wheeler, 
will achieve this (putting aside spurious, growing gauge 
modes |37|). In this gauge, the deviation vector zf grows 
quadratically with time, according to Eq. (60), which 
causes the second-order held h^ v (x; Zq) to li kewi se grow 
quadratically with time, as we see in Eq. (27b). Now 
suppose we perform a gauge transformation generated 
by = z /. This removes the first-order deviation, set¬ 
ting z'i = 0 and thereby eliminating the growing term in 
h'^(x-,Zo). But according to Eq. (121a), the transforma¬ 
tion generically introduces growing terms into h'^ v {x ; Zo). 
That is, the object’s motion away from Zq induces secu¬ 
lar growth in the second-order held; if we force the object 
to stay on z(( via a gauge choice, then we instead induce 
secular growth in the first-order held. 

This distinction between gauge and motion carries over 
to the self-consistent case, but the situation becomes 
more subtle. Here all of the helds h are function¬ 
als of a gauge-dependent worldline. Even for manifestly 
physical quantities such as the gravitational waveform, 
we hence seem to lack invariance. If we calculate, say, 
the Weyl scalar ipo from the linearized held in a partic¬ 
ular gauge,the result is invariant under the transforma¬ 
tion hr —?• but it is not invariant under 

fyU*] Kv\ z '] + 


To resolve this, we may hrst consider behavior of the 
approximation on short, order-unity timescales ~ e . 
On these timescales, the difference between /*! [z] and 
h^ v [z °] = h 1 is of second order. Hence, the hrst-order 
information in i/jq, as constructed from /i^[z], is invari¬ 
ant; there is merely an irrelevantly small bit of gauge- 
dependence buried within it. Similarly, if we construct 
a curvature invariant from /i* [*] plus /i^fz], it will be 
invariant through second order and contain a negligible 
third-order gauge-dependence. 

However, this description ignores the principal aim of 
the self-consistent expansion, which is to model long-term 
evolution. Specifically, look at an EMRI. A single orbit 
occurs on the timescale ~ e°, but the inspiral of the orbit 
occurs on a much longer timescale. Gravitational waves 
carry away orbital energy at the rate E/E ~ e, implying 
that the inspiral occurs on the time scale ~ E/E ~ 1/e. 
Hence, when modeling an EMRI we do not want to work 
on a fixed, e-independent interval of time; instead, we 
want to look at the limit e -> 0 in a domain V of size 
~ 1/e that blows up in the limit. 

Working on an e-dependent domain forces us to revise 
our thinking about gauge and motion. For our approx¬ 
imation to be sensible, we must require that on the do¬ 
main we work in, all excluded terms are much smaller 
than the included ones, and all included ones are well 
ordered. For example, to claim first-order accuracy, we 
must have z) ~ e° on T>, and we must have that all 


excluded terms are <C e°. These criteria strongly restrict 
the allowed choices of gauge. If, for example, we adopt a 
gauge in which there is no self-force, then /il„(ir; z) grows 
large with time; this is not an allowed gauge. Although 
these issues have not been studied in any detail, if we 
wish to stay within a class of gauges satisfying the above 
criteria, a natural requirement is to insist that all gauge 
vectors are uniformly small (and uniformly well ordered, 
such that e n+1 C+i < e n ££) in V. 

This implies that in well-behaved approximations in (D, 
the effects of the self-force are not pure gauge. The ac¬ 
celerated worldline z M eventually deviates from any given 
geodesic Zq by a very large amount, while with an al¬ 
lowable gauge transformation we may shift it only by a 
very small amount. In other words, although the self- 
forced deviation from Zq is pure gauge on a domain of 
size ~ e°, it is not pure gauge in a domain of size ~ 1/e. 
Furthermore, since the object is of size e and the gauge 
transformation may only shift the object by an amount 
of order e, within the allowed class of gauges, the gauge 
ambiguity in the object’s position is of the same order as 
the size of the object itself. For the purposes of model¬ 
ing long-term inspirals, this should not concern us; what 
concerns us is only the large, long-term changes in the 
object’s position. 

Let me summarize the practical relevance of all this. 
Suppose that in some gauge, one computes a metric per¬ 
turbation over a long timescale in the self-consistent ap¬ 
proximation, using the self-accelerated orbit as a source. 
Further suppose that one ends up with a well-behaved ap¬ 
proximation, in the sense that e n+1 h//f 1 <C e n /i" -C g^ v 
for all included terms, and that error estimates suggest 
that all excluded terms are smaller. Then one has in¬ 
cluded the correct invariant information about the orbit’s 
long-term evolution. One has also included some gauge- 
dependent information about the orbit, but on all scales, 
this information is negligibly small compared to the order 
of accuracy one should expect from one’s approximation. 
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Appendix A: Reparametrization transformations 

Section IIII Al describes how the coefficients in the ex¬ 
pansion 


z»{s,e) = Y J z n ^{s) (Al) 

n> 0 
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depend on the choice of parameter s. This appendix 
works out explicitly how the terms z£(s) transform under 
a reparametrization s —> s'(s, e). 

Let z' M (s',e) = z^(s(s', e), e) be the coordinates of the 
worldline in the new parametrization. Their expansion 
reads 

= ( A2 ) 
n> 0 

with coefficients z'£(s') := 9 (s', 0 )- 

We find the relationship between z and z^ by writing 


Appendix B: Identities on S 


In Sec. Ill A 2 we require various derivatives of quan¬ 
tities on a surface S = {z M (s, e)|s £ IR,e £ [0,oo)}. This 
appendix displays useful identities for that purpose. 

The surface can be generated by the two tangent vector 


fields v^ and zA = An important relation 

between these fields follows from the equality CiV^ = 
7^57 = since this implies 

= 0 = and hence 


(Bl) 


dx^(s',e) _ dz fl (s(s', e), e) 
de de 

dz lx {s,e) ^ ds(s',e) dz^(s,e) 


de de 

Evaluating at e = 0, we have 

ds 

L | -U 

+2 »s 


da 




e—0 


(A3a) 

(A3b) 


(A4) 


where Zq := From this we immediately find 

that the perpendicular piece of zf, unlike zf itself, is 
reparametrization invariant: 


g'M 

nx 


= z‘ 


ix- 


(A5) 


Similarly, we find the relationship between z% and z.f’ 
by writing 


d 2 z^(s',e) 

d? 


d 2 z^s(s',e),e) 

de 2 

<9 2 z M (s, e) d 2 s(s ', e) dz fl (s, 

de 2 de 2 ds 

0 ds(s\e) d 2 z^(s 1 e) 
de deds 

/ ds(s', e)\ 2 d 2 z^s,e) 

\ de ) ds 2 


(A 6 a) 

e) 


(A 6 b) 


Evaluating this at e = 0 and specializing to a coordinate 
system that is normal along 70 , we find 

U//J _ ^ 19 2 s dzg dsDz^ l(ds\ 2 D 2 zg 
‘ 2 N Zw + 2de 2 ds de ds 2\de) ds 2 ’ 

(A7) 


where the partial derivatives are evaluated at e = 0 . 
By projecting orthogonally to ito^, we find that unlike 
at first order, the perpendicular deviation Z 2 N j_ is not 
parametrization independent: 


It immediately gives us 


Using the identity (|Bl|), we can derive a host of others. 

(B2) 


ds (ds\ 


r _= 

V 


-) v a .pz a z^, 


dr \drJ 

since = 1/ yj—g^i^z". This in turn gives us 
(ds^\ 

C V K = — (z a Z^V/3- a + Z a V a — 2 KZ a V c 

+z“z /3 i 7 t) a;/ 3 7 ) , 
where k = —z a d n In 


Similarly, 


dr ’ 


C, v xA = 


Dv v 


= P» v 


where u" = ^ = f^z", and 
v 


dr ’ 
Dv v 


dr 


V„z“ = (i a ;/ 3 Z /3 ) ; ^u 7 

= v a + R a ^vH\ 


(B3) 

(B4) 

(B5) 


(B 6 a) 

(B 6 b) 


where the second line follows from Eq. (Bl) and the Ricci 
identity. 


Appendix C: Identities for gauge transformations of 
curvature tensors 

Let A[g\ be a tensor of any rank constructed from a 
metric g. (To streamline the presentation, I adopt index- 
free notation throughout most of this appendix.) Now 
define 


~ 2 NX — 2 2 NX 


ds Dz f 
de ds 


(AS) 


D 2 dz^ 

where I have used the geodesic equation —— — 


dz! 


= K 


ds 


However, after invoking Eq. (51) and Eq. (A4), we find 
that the quantity z^ is parametrization independent: 


*2f ~ z 2f 


(A9) 


6 n A[f u ...,f n ] := 

1 d n 

n\d\r • • • d\ n A[9 + Al/l + " ' + Kfn] k=o> ( C1 ) 

where “A= 0” stands for evaluation at A, = • • • = A n = 
0. The tensor-valued functional 6 n A[fi 1 ...,/„] is linear 
in each of its arguments /i, • • •, /„; it is also symmetric 
in them. In the case that all the arguments are the same, 
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we have S n A[h,... ,h) = ^ A[g + Xh] | A the piece 
of A[g + h] containing precisely n factors of h and its 
derivatives. 

The following identities are easily proved (and will 
be proved below) by writing Lie derivatives as ordinary 
derivatives: 


£ i A[g]=SA{^g\, (C2) 

\C\A\g] = \SA[C\g] + 8 2 A[^g, C^g], (C3) 

C e SA[h] = 8A[C^h] + 2 S 2 A[£ ig , h], (C4) 

where we can also write 8 2 A[£^g,h] = S 2 A[h, C^g] = 
\ ( S 2 A[h , C^g] + S 2 A[C^g, h ]). As an example, if A is the 
Ricci tensor, then 


CtRpvlg] = SR^g], (C5) 

\C\RM = \5R^[C 2 g]+8 2 R llv [C i g,Ci9\, (C6) 

Cc5R^[h\ = SR^fah] + 2 S 2 R^[h, C.g], (C7) 


where I have restored indices to avoid confusion with the 
Ricci scalar (though the same equation, with indices re¬ 
moved, would apply to the Ricci scalar). 

Before moving to the proof of these results, I note one 
of their consequences: when examining perturbations of 
a curvature tensor, we can derive transformation laws in 


two equally natural ways: directly from Eq. (881 or from 


the transformations of the metric perturbations. For ex¬ 
ample, in a vacuum background, Eq. (88) directly implies 


ASR^lh 1 ] = C il R lxv [g] =0, (C8) 

A [8R^[h 2 ] + S 2 R ltv [h 1 ,h 1 ]) = C^R^h 1 ]-, (C9) 


or the same equations can be found by instead applying 
Eq. ( [88] ) to the metric itself, writing 

A(SR^[h 2 ] + 8 2 R^[h\h 1 ]) 

= 8R^[h ,2 }+8 2 R^[h ,1 ,h n ] 

- {SR tlv [h 2 ] + 5 2 R„ v [h 1 ,h 1 ]) (ClOa) 

= 8Rfj. u [h~ +C^ 2 g+^C 2 1 g+C^ 1 h 1 } 

+ 8 2 R^ iy [h 1 +£^ 1 g,h 1 +C^ 1 g] (ClOb) 

= 8Rfw [£$2 d\ + ^8R tll j[C^ 1 g\ + SR^lC^h 1 ] 

+ 2 S 2 R^[h\C il g\+S 2 R^[C il g,C il g\ (ClOc) 


and then applying Eqs. (|C5j)-(|C7j) . 

Now, let us return to the proofs. To establish Eq. (C21, 
one can write the metric as a function of a parameter A 
along the flow generated by £ and then perform a Taylor 
expansion: 


Ci A \s\ = Tx A 


’<°> +A ! 


A—0 


= SA 


dg_\ 

RV A=0 


A=0 


= 5A[Ccg\. 


(Clla) 

(Cllb) 

(Cllc) 


Similarly, to establish Eq. (C3), one can write 




A=0 


C l A \-9] l dX 2 A 

= \8A[C 2 ^g\ + 8-AlC^g, C^g\, 


1 1 \ 2 dg 

+ 0 d\ 2 


A=0 


A—0 


(C12a) 


and to establish Eq. ( |C4[ ), one can write g as a function 
of parameters (A, e) along commuting flows, where h := 
I and then write 

de I e=0 ’ 




dXde 


-A 


g{X, 0 ) + e^(A, 0 ) 
dg , 


(C13a) 


A—e—0 

dg. 


9 ( 0 , 0 ) + A dX( 0 , 0) + e ^(°,0) 


+x X <o - a) 


A=e=0 


= 8A[£ ( h] + 2S 2 A[h,£ i g]. 


(C13b) 

(C13c) 


Appendix D: Gauge transformation of the 
second-order stress-energy 


According to Eq. (103), the gauge transformation of 
the (trace-reversed) second-order stress-energy, A T 2 


is 


given by Lie derivatives of the first-order stress-energy. 
This appendix evaluates those derivatives. 


The calculation can be performed following Sec. Ill B 2 


For any vector £(*, the Lie derivative T( tv (x; z) is given 
by Eq. |77| ) with the replacements and Up —> 

vA. The ordinary Lie derivative T[‘ 1 ' can be evaluated 
by following very similar steps as were used to derive 
Eq. j77]|, leading to 




a. ft 


= —m 


J&g$,{2u<- a, ^6(x,z) 

+ U a 'X S(x, z) 

- u a tf ± g^V 7 8(x,z)JdT, (Dl) 


where := a /£“ and £iy := . The total of 

the two Lie derivatives yields the simple result 


{£( + £ $) T r = ~ m J 9^9u'X 

' d i i|| 


dr 


+ £1 -p’) 8{x,z)dr. (D2) 


The ±^iV^i(s, z) terms in Eqs. (Dl) and ( |77| ) signal 
that the mass m is displaced from z M by an amount ; 
the lack of any such term in Eq. (D2) signals that the 


displacements due to the two derivatives precisely cancel 
one another, leaving the mass m moving on . 
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Equation (103) involves derivatives of T* , not of Tf 1 '. 


So we have 


£$,i)Ti + 2(C( m;q )5^/3 
+ — 2^(p.',v)9ot0 ~ 29^^ a -p)Ti ^) (D3) 


with (£^ x + £^)T^ given by Eq. (D2). 

These are the results for the self-consistent case. For 
the Gralla-Wald case, we simply drop the term and 


set z^ = Zq: 


— {g^agup 29^’ J 9ap)£-^ 1 Ti + 2(£( /[i . a )<?i//3 


' 9p-at,( v -p) 2^(n;v)9a.0 \9 pvZ, 0 c,fi)'£'l (D4) 


where £ f l T"' 3 is given by Eq. ( |PT| ) with z^ zft. Un¬ 
like in the self-consistent case, where the V M 6(a;, 2 ) terms 
cancelled in the final result, in the Gralla-Wald case there 
is a term oc 2 ), corresponding to the center of 


mass having been displaced by an amount Az^ 


relative to z[ 
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